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Abstract. We study three families of labelled plane trees. In all these trees, the root is 
labelled 0, and the labels of two adjacent nodes differ by 0, 1 or —1. 

One part of the paper is devoted to enumerative results. For each family, and for all 
j G N, we obtain closed form expressions for the following three generating functions: 
the generating function of trees having no label larger than j; the (bivariate) generating 
function of trees, counted by the number of edges and the number of nodes labelled j; 
and finally the (bivariate) generating function of trees, counted by the number of edges 
and the number of nodes labelled at least j. Strangely enough, all these series turn out 
to be algebraic, but we have no combinatorial intuition for this algebraicity. 

The other part of the paper is devoted to deriving limit laws from these enumerative 
results. In each of our families of trees, we endow the trees of size n with the uniform 
distribution, and study the following random variables: M„, the largest label occurring in 
a (random) tree; Xn{j), the number of nodes labelled j; and X^{j), the number of nodes 
labelled j or more. We obtain limit laws for scaled versions of these random variables. 

Finally, we translate the above limit results into statements dealing with the integrated 
superBrownian excursion (ISE). In particular, we describe the law of the supremum of 
its support (thus recovering some earlier results obtained by Delmas), and the law of its 
distribution function at a given point. We also conjecture the law of its density (at a 
given point). 



1. Introduction 

We study in this paper three famihes of labelled plane trees. In all these trees, the root 
is labelled 0, and the labels of two adjacent nodes differ by 0, 1 or —1. 

More precisely, the first family we consider is the set of plane trees, and the increments 
of the labels along edges are constrained to be ±1. In the closely related second family, 
these increments can be 0,±1. The third family is a bit different. It is simply the set of 
(incomplete) binary trees, in which the nodes are labelled in a deterministic way: the label 
of a node is the difference between the number of right steps and the number of left steps 
occurring in the path that yields from the root to the node under consideration. See Figure^] 
for an illustration. We call this labeUing the natural labelling of the binary tree. Note that 
the label of each node is simply its abscissa, if we draw the tree in the plane in such a way 
the right (resp. left) son of a node lies one unit to the right (resp. left) of its father. For this 
reason, we will sometimes call these labelled binary trees naturally embedded binary trees. 
More generally, for any plane labelled tree, we may consider that the label of each node tells 
where to embed it in Z; hence the title of the paper. 

In each of these three famihes, we endow the set of trees having a given size (say, n edges) 
with the uniform distribution. We address (via generating functions) the following three 
questions: 
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Figure 1. A labelled plane tree with increments ±1. — A labelled tree 
with increments 0,±1. — A naturally embedded binary tree. 

(1) What is the maximal label that occurs in the tree? This label is in fact a random 
variable M„. We prove that Mnf-n}/^ converges in distribution to a random variable 
N having a density. We give this density explicitly. We also compute the moments 
of N and prove the convergence of the moments of Mn/n}^'^ to those of N . 

(2) How many nodes of the tree have label j? Let X„(j) denote the corresponding 
random variable. If j is fixed, and n goes to infinity, then the answer to this question 
is independent of j. We prove that for any j € Z, the variable Xn{j)/n^^^ converges 
in distribution to cT~^/^, where c is a constant depending on which family of trees 
we consider, and T follows a unilateral stable law of parameter 2/3. 

Given that the maximal label grows Hke n^/"*, we get a better insight on the 
label distribution by asking how many nodes in a tree of size n have label [An^/^J . 
We prove that, for any A S R, the random variable Ar„([An^/"'J converges 
in distribution to a limit variable ^(A). This variable admits a Laplace transform, 
which we give explicitly. The convergence of the Laplace transform, and of the 
moments, hold as well. We say we have obtained a local limit law for embedded 
trees, because we look at one value of the labels only. 

(3) Finally, we also obtain a global limit law by studying the variable X^{j) that gives 
the number of nodes having label j at least. Remarkably, we prove that X+(0)/n, the 
(normalized) number of nodes having a non-negative label, converges to the uniform 
distribution on [0, 1]. More generally, for A £ M, the variable X^{Xn^^'^)/n converges 
in distribution to a variable Y~^{X). This variable admits a Laplace transform, which 
we give expHcitly. Once again, the convergence of the Laplace transform, and of the 
moments, hold as well. 

The laws of A^, Y{X) and Y~^{X) naturally depend on which family of trees we consider, but 
only by a simple normalization factor. 

1.1. Embedded trees and the integrated superBrownian excursion 

Why should one study such labelled trees? 

The first two classes of trees we consider have a close connection with certain families of 
planar maps jHllHlE]- In particular, the diameter of a random quadrangulation having n 
faces is distributed like the largest label in non-negative random trees of our second family. 
Moreover, once scaled by n^/'*, this diameter has the same limit law as (M„ — m„)n~^/'*, 
where M„ (resp. m„) is the largest (resp. smallest) label occurring in a random tree of our 
second family 

The third class we study is the good old family of binary trees, and this may suffice to 
motivate its study! More seriously, the three questions addressed above have, for binary 
trees, a natural geometric formulation. The random variable M„ (the maximum label) tefis 
us about the "true width" of a binary tree (as opposed to the maximal number of nodes lying 
at the same level, which is known to grow like ^/rl). More generally, the variables Ar„(j) tell 
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Figure 2. An (incomplete) binary tree having horizontal profile [1, 2, 4, 3, 2] 
and vertical profile [2, 2, 4; 2, 1, 1]. 



us about the vertical profile of the tree (as opposed to the horizontal profile which describes 
the repartition of nodes by level 53)- See Figure IS 

We may also invoke an a posteriori justification to the study of these trees: the form of 
the generating functions we obtain is remarkable, whatever family of trees we consider, and 
suggests that there must be some beautiful hidden combinatorics in these problems, which 
should be explored further. 

However, the main motivation for this work is the connection between embedded trees 
and the integrated superBrownian excursion (ISE). Choose one of the three families of trees, 
and consider the following random probability distribution on M: 

M'l = ^n(j)^cjn-i/4, (1) 

n + 1 ^-^ ■' 

where Xnij) is the (random) number of nodes labelled j, 5x denotes the Dirac measure at x, 
and the constant c equals -\/2 for the first family, -\/3 for the second one and 1 for the family 
of binary trees. Then ^„ is known to converge weakly to a limiting random probability 
distribution called the ISE [HEBIESEni- See Figure for simulations of /^„. 

Our limit results provide some information about the law of the ISE. For instance, we 
prove that cM„n~^/^, the largest point having a positive weight under converges in law 




Figure 3. The plot of Xn{i) vs. j for random binary trees with n = 1000 nodes. 
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to NisE, the supremum of the support of the ISE. We denote this by 

The results we obtain for the limit law of Af„n~^/^ thus translate into expressions of the 
moments, distribution function and density of the supremum of the ISE. Note that the 
moments were already obtained by Delmas 22j- Our second limit result deals with the 
random variables X„([An^/'*J ). Observe that 

/i„(cA-cn-i/4^cA] = ^X„(LAni/4j). (2) 
n + 1 

This leads us to conjecture that the random variable Y{X) involved in our local limit law 
satisfies 

Y{X) - c/ise(cA) (3) 
where /ise is the (random) density of the ISE. Similarly, 

/i„[cA,+cx)) = ^X+(rAni/4]), 
n + 1 

and we prove that the random variable ^^(A) involved in our global limit law satisfies 

where gisE is the (random) tail distribution function of the ISE. The results we obtain about 
the laws of F (A) and Y'^ (A) thus translate into formulas for the Laplace transforms of /ise(A) 
and (7ise(A) (the formula for /ise (A) being conjectural). 

Our conjecture on /ise is naturally supported by the fact that the law of Y{X/c)/c is 
independent of the tree family we start from. This is one of the reasons why we consider as 
many as three families of trees. The other reasons involve the connections with planar maps, 
the remarkable form of the generating functions we obtain, and our unshakeable interest in 
binary trees. The details of the calculations are only given for the first of the three families 
(Sections 121 to while the results are merely stated for the other two families (Section EJ. 

Let us finally mention that the moments of the center of mass of the ISE have recently 
been determined by two different approaches [3 El ■ In our discrete setting, this boils down 
to studying the convergence of the variable 

1.2. Overview of the paper 

The starting point of our approach is a series of exact enumerative results dealing with 
our first class of trees: plane trees in which the labels of adjacent nodes differ by ±1. These 
results are gathered in the next section. We obtain for instance an expHcit expression for 
the bivariate generating function of labelled trees, counted by the number of edges and the 
number of nodes labelled j (for j fixed). This section includes, and owes a lot to, some 
results recently obtained by Bouttier, Di Francesco and Quitter 01^1 on the enumeration of 
trees having no label greater than j. This part of our work raises a number of challenging 
combinatorial questions — why are these expressions so simple? — which are not addressed 
in this paper. 

The limit behaviours of the random variables Af„, X„([An^/'*J) and X^{Xn^^'^) are re- 
spectively established in the next three sections (Sections to |^. The main technique that 
we use is the "analysis of singularities" of Flajolet and Odlyzko It permits to extract 
the asymptotic behaviour of the coefficients of a generating function. This technique has 
already proved useful in numerous occasions, in particular for proving Hmit theorems that 
are similar in fiavour to the ones obtained in this paper: these theorems deal with the height 
of simply generated trees and their profile, which are known to be related to the height of 
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the Brownian excursion and its local time \W\ 113^ . This technique is carefully exemplified 
in Section El (which is devoted to the maximal label) before the more difficult questions of 
the local and global limit laws are attacked (Sections ^ and 1^. 

Finally, two other families of trees are briefly studied in Section El trees with increments 
0,±1 and naturally embedded binary trees. The emphasis is put on their enumerative 
properties, which turn out to be as remarkable and surprising as those of our flrst family of 
trees. The limit laws we obtain are (up to a scalar) the same as for the first family. 

Let us conclude with some notation and a few definitions on formal power series and 
generating functions. Let K be a field. We denote by K[t] the ring of polynomials in t with 
coefficients in K, and by K(t) the field of rational functions in t with coefficients in K. We 
denote by IK[[i]] the ring of formal power series in t with coefficients in K. If A{t) G lK[[t]] 
and n e N, the notation [t"]A(t) stands for the coefficient of in A{t). The series A{t) 
is said to be algebraic over K(t) if it satisfies a non-trivial polynomial equation of the form 
P{t,A{t)) — 0, where P is a bivariate polynomial with coefficients in IK. In this case, the 
degree of A{t) is the smallest possible degree of P (in its second variable). 

Let ^ be a set of discrete objects, equipped with a size that takes nonnegative integer 
values. Assume that for all n g N, the number of objects of A of size n is finite, and denote 
this number by a„. The generating function of the objects of A, counted by their size, is the 
formal power series 

A{t) = 

n>0 

The above notions generalize in a straightforward way to multivariate power series. Such 
series arise naturally when enumerating objects according to several parameters. 

2. Enumerative results 

We consider in this section (and in the three following ones) our first family of labelled 
plane trees: the root is labelled 0, and the labels of two adjacent nodes differ by ±1. 

2.1. Trees with small labels 

The first enumerative problem we address has already been studied by Bouttier, Di 
Francesco and Guitter OEj. It deals with the largest label occurring in a tree. For j e N, 
let Tj = Tj{t) be the generating function of labelled trees in which all labels are less than 
or equal to j. The indeterminate t keeps track of the number of edges. Let T = T{t) be the 
generating function of all labelled trees. Clearly, Tj converges to T (in the space of formal 
power series in t) as j goes to infinity. It is very easy to describe an infinite set of equations 
that completely defines the collection of series Tj . 

Lemma 1. The series T satisfies 

T = 1 + 2tT'^. (4) 

More generally, for j > 0, 

= i + t(r,_i + r,+i)T, 

while Tj = for j < 0- 

Proof. The two ingredients of the proof will be useful for the other enumerative problems 
we address below. Firstly, replacing each label khy j — k shows that Tj is also the generating 
function of trees rooted at j and having only non-negative labels (we say that a tree is rooted 
at j if its root has label j). Secondly, consider such a tree and assume it is not reduced to a 
single node. The root has a leftmost child, which is the root of a labelled subtree, rooted at 
j ± 1 and having only non-negative labels. Deleting this subtree leaves a smaller tree rooted 
at j, having only non-negative labels (see Figure 0J. The result follows. 
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Figure 4. The decomposition of plane labelled trees. 



The above lemma shows that the series T, counting labelled trees by edges, is algebraic, 
and the short proof we have given provides a simple combinatorial explanation for this 
property. What is far less clear — but nevertheless true — is that each of the series Tj is 
algebraic too, as stated in the proposition below, which we borrow from |^|^. These series 
will be expressed in terms of the series T = T{t) and of the unique formal power series 
Z = Z{t), with constant term 0, satisfying 

Observe that T and Z are related by: 

(6) 

Proposition 2 (Trees with small labels OE]). LetTj = Tj{t) he the generating function 
of trees having no label greater than j. Then Tj is algebraic of degree (at most) 2. In 
particular, 

Moreover, for all j > —1, 



To = 1 - lit - + 4i (3 + 2t)ro - 16<^To^ 

T = T - — - (7^ 

^ (l-Zi+2)(l-ZJ+4)' ^ ' 

where Z = Z{t) is given by 

Proof. It is very easy to check, using that the above values of Tj satisfy the recurrence 

relation of Lemma and the initial condition T^i = 0. How to discover such a formula 
is another story, which is told in 0. The remarkable product form of Tj still awaits a 
combinatorial explanation. 

The equation satisfied by Tq is obtained by eliminating T and Z from the case j — 
of 0. Then an induction of j, based on Lemma impHes that each Tj is quadratic (at 
most) over ' ' 



Remarks 

1. The product form 10, combined with the facts that T is quadratic over Q(t) and Z is 
quadratic over Q(r), shows that Tj belongs to an extension of Q(<) of degree 4. This is true, 
but not optimal, since Tj is actually quadratic over Q(i). Hence this product form does not 
give the best possible information on the degree of Tj . 

2. The trees counted by Tq (equivalently, the trees having only non-negative labels) are 
known to be in bijection with certain planar maps called Eulerian triangulations [Ol- Through 
this bijection, the number of edges of the tree is sent to the number of black faces of the 
triangulation. These triangulations are nothing but the dual maps of the bicubic (that is. 
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bipartite and trivalent) maps, which were first enumerated by Tutte jJE]- In particular, the 
coefficients of ro(i) are remarkably simple: 

n>l / \ / 

2.2. The number of nodes labelled j 

Let us now turn our attention to a bivariate counting problem. For j e Z, let Sj = Sj{t, u) 
be the generating function of labelled trees, counted by the number of edges (variable t) and 
the number of nodes labelled j (variable u). Clearly, Sj{t, 1) = T{t) for all j. Moreover, an 
obvious symmetry entails that Sj ~ S-j. 

Lemma 3. For j ^ 0, 

S, = l + t{Sj-i+S,+i)S, (8) 

while for j = 0, 

So=u + t{S-i + Si)So=u + 2tSiSo. (9) 

Proof. Observe that Sj = Sj{t, u) is also the generating function of labelled trees rooted at 
j, counted by the number of edges and the number of nodes labelled 0. The decomposition 
of trees illustrated in Figure then provides the lemma. The only difference between the 
cases j = and j lies in the generating function of the tree reduced to a single node. 

■ 

Again, the series Sj{t,u) turn out to be algebraic, for reasons that currently remain 
mysterious (from the combinatorics viewpoint). They can be expressed in terms of the series 
T and Z given by (|3El- 



Proposition 4 (The number of nodes labelled j). For any j G Z, the 

function Sj = Sj{t,u) that counts labelled trees by the number of edges and the number of 
nodes labelled j is algebraic over Q(T, w) of degree at most 3 (and hence has degree at most 
6 over Q{t, u) ). More precisely, 

{T-S^ ^ J2^-T-^ (10) 
and all the Sj belong to Q{t,u, Sq). Moreover, for all j > 0, 

^ ji+M£^)a+M^^ 

' (l + /iZ^+i)(l+/iZJ+3)' ^ ) 

where Z = Z{t) is given by ^ and ^ = /i(i, u) is the unique formal power series in t 
satisfying 

^ (l + Z^)(l + /iZ)(l + /iZ^)(l + /iZ3) 

^ ^ '{l + Z){l + Z + Z^){l-Z)^{l^fiZ'^)- ^ ' 

The series ^i{t,u) has polynomial coefficients in u, and satisfies ^{t, 1) = 0. It has degree 3 
overQ{Z,u) and 12 overQ{t,u). 

At some point, we will need a closed form expression for /i in terms of Z. Here is one. 

Proposition 5. Write 

(m-1)Z(1 + Z2) 

V — 



{l + Z){l + Z + Z^){l- Zf 
Then the algebraic series fi involved in the expression Hll|i of Sj, and defined by Ijl2|l . is 

, , 1 / 2 \ 

\i{t,u) = 



Z"^ I 1 + i;(l - Z)2/3 + 2/3^/3 + f2(l - zy cos((/)/3) 
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where 

-9v{l + AZ + Z^) + v^{l - ZY 



(3 + i;2(l-Z)4)3/2 

Proof of Propositions ^ and [sj^ First, observe that the family of series Sq, Si, S2, ... is 

completely determined by (|Hl (taken for j > 0) and the second part of The fact that 
for any series /x £ Q(w)[[t]], the expression ifTTll satisfies © for all j > is a straighforward 
verification, once t and T have been expressed in terms of Z (see ^ and ©). The form 
of Ijllll is borrowed from 0. In order for to be the correct expression of Sj, it remains 
to satisfy the second part of Q . This last condition provides a polynomial equation relating 
fj., T, Z, t and u. In this equation, replace t and T by their expressions in terms of Z (given 
by (EHEl)- This gives exactly l(T2|l . It can be easily checked that /.i has degree 6 over Q(T, u) 
and degree 12 over Q(i, u). 

The equation satisfied by Sq is obtained by eliminating /i and Z (using Ijl2|l and 
from the expression lllip of Sq. This equation gives an equation of degree 6 over Q(t,w) if 
one eliminates T thanks to 

Now the equations ©, © and Q, combined with an induction on j, imply that for j > 1, 
the series Sj belongs to the field Q(r, u, S'o), which has just been proved to be an extension 
of Q(r, u) of degree 3. This concludes the proof of Proposition^! 

Let us finally prove Proposition The equation Ijl2|l that defines jjL can be rewritten 

_ V jl + fiZ){l + fiZ^){l + nZ^) 
Z 1 - ^iZ^ 

Hence /i is the unique formal power series in v (with rational coefficients in Z) that satisfies 
the above equation and equals when u is 0. It is not hard to check that the closed form 
expression we give satisfies these two conditions. 

■ 

Remarks 

1. The product form ((TTll of Proposition^lrefines the product form (0 that deals with trees 
with small labels. Indeed, when u = 0, Eq. l[T2ll gives fj. = —1, and the expression of Sj{t, 0) 
coincides, as it should, with the expression of T,_i(i) given by Proposition El 

2. There exists an alternative way to derive an equation for Sq from the system of Lemma^l 
As was observed in [£| p. 645] for the problem of counting trees with bounded labels, Eq. © 
implies that for .7 > 1, 

I{Sj-i,Sj) = I{Sj,Sj+i) 
where the "invariant" function / is given by 

I{x, y) = xy{l - tx){l - ty) + txy -x-y. 
But Sj converges to T as j goes to infinity, in the set of formal power series in t. This implies 

I{Sa,Si)=I{T,T). 

Eliminating Si between the above equation and Q gives an equation between S'o, T and t. 

2.3. The number of nodes labelled j or more 

Let us finally study our third and last enumeration problem. For j £ Z, let Rj = Rj{t, u) 
be the generating function of labelled trees, counted by the number of edges (variable t) and 
the number of nodes labelled j at least (variable u) . 



^AU the calculations in this paper have been done using Maple. We do not recommend the reader to 
check them by hand. 
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Lemma 6. The set of series Rq, Ri, R2, . ■ . is completely determined by the following equa- 
tions: for j > 1, 

Rj = l + tRj + R^+i) (13) 

and 

Ra{t,u) =uRi{tu,l/u). (14) 

More generally, for all j G one has: 

R^j(t,u) ^ uRj^i{tu,l/u). (15) 

Proof. For all j G Z, the series Rj = Rj{t, u) is also the generating function of trees rooted 
at j, counted by their number of edges and the number of nodes having a non-positive label. 
The equation satisfied by j, for j > 1, follows once again from the decomposition of trees 
illustrated in Figure ^ It remains to prove the symmetry relation (I15|l . For any tree r, 
let n<o(T) denote the number of nodes of r having a non-positive label. We use similar 
notations for the number of nodes having label at most j, etc. Let 7^^„ denote the set of 
trees rooted at j and having n edges. As observed above, 



n>0 TeT-j^n n>a TeT_j,„ 



because a tree with n edges has a total of ?i + 1 nodes. A translation of all labels by —1 gives 

R.,{t,u)^uY,{tur Yl """-"'"^ 

n>Q TeT_j_i.„ 

while replacing each label k by —k finally gives 

R-jit,u) = u^(t?i)" M""so(r) ^ uRj+i{tu, 1/u). 

n>0 TeTj + i,„ 



Again, the series Rj are algebraic, and admit a closed form expression in terms of T and 

Z. 

Proposition 7 (The number of nodes labelled j or more). Let j £ Z. The generating 
function Rj (t, u) = Rj that counts labelled trees by the number of edges and the number of 
nodes labelled j or more is algebraic of degree at most 2 over Q{T (t) , T (tu)) . Hence it has 
degree at most 8 over Q{t,u). More precisely, it belongs to the extension of Q{T(t),T{tu)) 
by 



(T + f )2 - ATf{T - l)(r - 1) 



where T = T{t) and T — T{tu). 
Moreover, for all j > 0, 

{l + iyZ^){l + i'Zi+^) 
(l + :/ZJ'+i)(l + ^ZJ'+3)' 

where Z = Z(t) is given by Q and v = v{t, u) is a formal power series in t, with polynomial 
coefficients in u, which is algebraic of degree 4 over Q(u, Z), and of degree 16 over Q{t, u). 
This series satisfies v{t, 1) = 0. The first terms in its expansion are: 

v{t, u) = (?i - 1) (^1 + 2 + (7 u + 6 u^) + (32 m + 36 + 23 u^) + 0(t^)) . 

Before we prove this proposition, let us give something like a closed form for v. Since v 
has degree 4 over Q(u, Z), and Z has degree 4 over Q(<), the series v is in theory expressible 
in terms of radicals... It turns that this expression is less terrible than one could fear. 
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Proposition 8. Define the following four formal power series in t with polynomial coeffi- 
cients in u: 

= o(t,u) = 1 — 8(u — 1)^ -;- = , 



V=V{t,u) 



V~5 



1 - 



l-8tM 

l-8t 



and 



A(i,u) = (l-F) 
P={l + Z) 



4 

AZV^ 

(X+zy- 



2VZ 

Then P has degree 16 over Q(t,u), degree 2 over Q{V,Z), and satisfies the following "La- 
grangian" equation: 

P=-I—{l+P){l + ZP). 
Moreover, the algebraic series v involved in the expression Ijl6|l of Rj is 

_ P 1 - P(l + Z) - p2(i + z + Z^) 
^ ^ Z 1 + Z + Z2 + PZ{1 + Z)- P^Z^ " 

Proof of Proposition 13 We have already checked, in the proof of Prepositional that for 
any formal power series in t, the series defined by ltT6|l for j > satisfy the recurrence 
relation f^l'M for j > 1. It remains to prove that one can choose v so as to satisfy H14II . 
For any formal power series ^ in i having rational coefficients in u, we denote by A the 
series A(t, u) — A{tu, 1/u). Observe that A = A. With this notation, if Rj is of the generic 
form Ijl6|l , the relation Ijl4|l holds if and only if 

^ ^. _ ^ r (1 + '^m + i.z^){i + z>z)(i + DZ') 

T (l + i/Z4)(l + i>Z2)(l + j>Z4) ■ ^ ' 

Let MmM denote the space of polynomials in u, with real coefficients, of degree at most 
m. Let E„[M][[t]] denote the set of formal power series in t with polynomial coefficients in u 
such that for all m < n, the coefficient of has degree at most m. Observe that this set of 
series in stable under the usual operations on series: sum, product, and quasi-inverse. Write 
u = J2n>o '^niu)t"'. We are going to prove, by induction on n, that II17|I determines uniquely 
each coefficient Un{u), and that this coefficient belongs to M„+i[m]. 

First, observe that for any formal power series the right-hand side of Ijl7|l is u + 0{t). 
This implies vq{u) ~ u — \. Now assume that our induction hypothesis holds for all m < n. 
Recall that Z is a multiple of t: this implies that i^Z belongs to K„[M][[t]]. The induction 
hypothesis also implies that the coefficient of in ui> belongs to Km+i[u], for all m < n. 
Note that Z = Z{tu) = tu-\-0{t^) is a multiple oft and u and also belongs to R„[u][[t]]. This 
implies that vZ belongs to R„[u][[t]] too. The same is true for all the other series occurring 
in the right-hand side of ifTTji . namely T, T, Z, Z. Given the closure properties of the set 
R„[u][[t]], we conclude that the right-hand side of ifTTji . divided by m, belongs to this set. 
Moreover, the fact that Z and Z are multiples of t guarantees that the coefficient of i" in 
this series only involves the Vi(u) for i < n. By extracting the coefficient of t" in IjlTII . we 
conclude that f„(M) is uniquely determined and belongs to C K„+i[m]. 

This completes the proof of the existence and uniqueness of the series f satisfying (tT7|l . 
Also, setting u = 1 (that is, T = T and Z = Z) in this equation shows that 1) = 0. 

Let us now replace t by tu and u by l/w in ifTTji . This gives: 

^ 1 ^ {l + i>Z){l + i>Z^)il + i^Z){l + ,.Z^) 

uf {l + uZ^){l + uZ^){l + uZ^) ■ ^ ' 
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In the above two equations, replace T by its expression ij^J in terms of Z . Similarly, replace T 
by its expression in terms of Z . Finally, it follows from Q and from the fact that Z ~ Z{tu) 
that 

.^^ii±£ni±l!l. (19) 

^ (1 + Z)4(1 + Z2) ^ ' 

Replace u by this expression in (I17|l and Ijl8|l . Eliminate v between the resulting two equa- 
tions: this gives a polynomial equation that relates v, Z and Z , of degree 2 in v. The 
elimination of Z between this quadratic equation and ifTfljl provides an equation of degree 
4 in f that relates i/ to Z and u. Finally, the elimination of Z shows that v is algebraic of 
degree 16 over Q{t,u). 

Let us now focus on the first part of the proposition. From the form Hl(i|l . and the fact 
that f has degree 4 over Q(u, Z) and Z has degree 4 over Q{t), we conclude that the degree 
of Rj over Q{t, u) is a divisor of 16. Let us prove that is is actually a divisor of 8. The proof 
goes as follows: 

(1) Using the generic form lfT6jl . and the equations satisfied by T, Z and i^, we obtain a 
polynomial equation of degree 8 over Q{t, u) for Rq. 

(2) Using Q to express t in terms of T, and 

_ T2 1 - f 

" ~ T2 1 — r' 

(which also follows from 1^), we convert the equation satisfied by i?o into a polyno- 
mial equation (still of degree 8 in Rq) relating Rq to T and T. This equation factors 
into four quadratic polynomials in Rq. The factor that actually vanishes is identified 
by setting u = 1 (in which case T = T = Rq). 

(3) From this equation, we conclude that i?o belongs to the extension of Q{T, T) gener- 
ated by 

VAi = \J{T + ff - ATf{T - l){f - 1). 

Observe that this extension of Q(t, u) is left invariant by the transformation A. 

(4) From the fact that Ri = uRq (see dJ), we conclude that Ri also belongs to 
Q(T,f,^/Al). 

(5) The recurrence relation ifT^jl on the Rj allows us to extends this to all Rj, for j > 0. 

(6) Finally, l|T5| shows that our algebraicity result actually holds for all Rj, for j G Z. 

■ 

Proof of Proposition |8j In the course of the proof of Proposition we have obtained 
a polynomial equation P(i/, Z, u) = 0, of degree 4 in relating the series i/(t, u),Z{t), and 
the variable u. This equation is not written in the paper (it is a bit too big), but it follows 
from Ijl7|l and H18II . In this equation, replace u by its expression in terms of S and Z. Then 
replace S by its expression in terms of V: the resulting equation factors into two terms! Each 
of them is quadratic in v. In order to decide which of these factors cancels, one uses the fact 
that when u = 1 (that is, V = 0), the series v must be 0. It remains to solve a quadratic 
equation in i/. Its discriminant is found to be A, and one may find convenient to introduce 
the series P which is Lagrangian in V. 

■ 

Remark. Again, the product form ((T6ll of Proposition includes as a special case the 
enumeration of trees with labels at most j — 1, obtained when u = 0. Indeed, (I17|l shows 
that u ^ —1 when u ~ 0, and (|16ll then reduces to 0. 
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3. The largest label, and the support of the ISE 

Let To denote the set of labelled trees (rooted at 0), and let 7o „ denote the subset of 
To formed by trees having n edges. We endow 7o,„ with the uniform distribution. In other 
words, any of its elements occurs with probability 

1 

where C„ = is the nth Catalan number, and is well-known to be the number of 

(unlabelled) plane trees with n edges. 

Let Mn denote the random variable equal to the largest label occurring in a random tree 
of 7o_„. The law of Af„ is related to the series Tj studied in Proposition |2l 

P(M„<j) = |^. 
Let us define a normalized version of A'/„ by 

The aim of this section is to prove the convergence of Nn in distribution^. 

Theorem 9. As n goes to infinity, the random variable Nn converges in distribution to a 
non-negative random variable N . The tail distribution function of N , defined by G(A) = 
F{N > A), satisfies 

4 

12 f v^e^ , 6 1 — cosucoshu = .^/ux*-! , 

G(A = / 5 dv = / ; e" 'du 

W'^ Jr sinh (Aw) V'^A*' Jo (cosh u - cosu)^ 

where the contour T is formed of two half-lines: 

r = {1 -^e-"^/^^ G (oo,o]}u{l + ^e-"^/^^ e [o,oo)}. 

Equivalently, the variable N has density 



f(X) = ^ f C0Sh(Az;)w6^" ^ p l-C0S^C0Shu ^3 4 _ -nV(4A*)rf^ 

iy^ Jr sinh3(Au) V^A" 7o (coshu - cosm)^ 

with respect to the Lebesgue measure on IR+. The moments of N are finite, and admit simple 
expressions: 

E(vv) = ^^(^, m')^3v^, 

and for fc > 3, 

24^k\C{k - 1) 



lU 



2'=r((A:-2)/4) • 

Finally, the moments of Nn = Mn/n^^'^ converge to the moments of N. 
The functions G and / are plotted in Figure 

The proof of this theorem will be split into four subsections fSections 13. II to . In view of 
the following proposition, this theorem gives the density, distribution function and moments 
of the supremum of the support of the ISE. 



■^The above convention will be used throughout the paper: if a random variable depending on n is 
denoted by some letter of the alphabet, then its suitably normalized version is denoted by the next letter of 
the alphabet. 
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Proposition 10 (The supremum of the support of the ISE). Let Nigs denote the 
supremum of the support of the ISE 

iVisE = sup{y : ^isE(y,oo) > 0}. 
Then iVisE has the same law as the random variable ^/2N described in Theorem^ 
Remark. The moments of A'^ise are thus 

3^/^ 
V2r(3/4)' 
and for fc > 3, 

V2 r((fc-2)/4) 

They were aheady obtained by Delmas jTJ] using a completely different (and continuous) 
approach. The expressions he gives actually differ from ours by a factor 2'^/^, due to a 
different choice of normalization. Note that the zeta function also appears in the moments 
of the maximum of the Brownian excursion, which follows a theta law ^OI- This law is known 
to describe the limiting normalized height of simple trees ^Hj. Finally, let us mention that 
another, more compHcated expression of the density of the limiting variable N was obtained 
in (maybe in a sHghtly less rigorous way). Proposition E| is proved in Section IT5l 

3.1. Convergence of the distribution function 

We prove in this section that the tail distribution function of 7V„ converges pointwise. Let 
A > and j = [An^/'*J . The probability we are interested in is 

P(7V„ > A) = P(A/„ > Ani/4) = P(M„ > j) - ^^^iP^' (^0) 

where 

_ (1 + Z)' ZJ+i (1 + Z + Z2) (1 - Z)' 
it) = U,-T- Tj - (1 ^ ^2) (1 „ ^,+2) (1 _ z^+i) (^^^ 

is the generating function of trees having at least one label greater than j. This algebraic 
series has a positive radius of convergence^, and by Cauchy's formula. 



2iTr Jc 

1 f {i + zfzi+^(i + z + z^){i-zf dt 



('221 

2i'K Jc (1 + ^2) (1 - (1 _ ' ^ > 

for any contour C included in the analyticity domain of Uj and enclosing positively the origin. 

This leads us to study the singularities of Uj, and therefore those of Z. We gather in the 
following lemma a few properties of this series. 

^So do all algebraic power series 
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Lemma 11 (Analytic properties of Z). Let Z = Z{t) he the unique formal power series 
in t with constant term satisfying Q. This series has non-negative integer coefficients. 
It has radius of convergence 1/8, and can he continued analytically on the domain V = 
C\ [1/8, +oo). In the neighborhood of t = 1/8, one has 

Z(i) = 1 - 2(1 -80^/^ + 0(\/T~~8i). (23) 

Moreover, \Z{t)\ < 1 on the domain V . More precisely, the only roots of unity that are 
accumulation points of the set Z{V) are 1 and —1, and they are only approached hy Z(t) 
when t tends to 1/8 and when \t\ tends to oo, respectively. 

Proof. In order to establish the first statement, we observe that 

Z = W{1 + Zf 

where W = W{t) is the only formal power series in t with constant term zero satisfying 

W = t + 2W^. (24) 

These equations imply that both W and Z have non-negative integer coefficients. 

The general approach for studying the singularities of algebraic series (see for instance |18) 1 
gives the second part of the lemma (up to lf23jl V The polynomial equation defining Z{t) has 
leading coefficient t and discriminant 4(1 — 8^)"^, so that the only possible singularity of Z is 
1/8. Alternatively, one can exploit the following closed form expression: 



^/T^ [ Vi - 4t + v^r^^ - ^/2(l - uyi'^ 



Z{t) = ^ (25) 

Let us now come to the third part of the lemma, and prove that \Z{t)\ never reaches 1 on 
the domain V. Assume Z{t) = e**, with 9 g [— 7r,7r]. Prom one has 

1 COS0 1 /, / N 

t^te where to ^ and t^e(-7r, tt). 

8 cos"*(w/2) 

This shows that t is real, and belongs to (— oo, 1/8). But the expression (l25|l of Z[t) shows 
that Z{t) is real, which contradicts the hypothesis Z{t) = e*^, unless = 0. But then i = 1/8 
and does not belong to the domain V. Hence the modulus of Z never reaches 1 on V. One 
can actually prove that, for 6 e (— 7r,0), 

cos u 

but we do not need so much precision here. 

Finally, if a sequence f„ of V is such that Z[tn) e'^ as n ^ oo, with 9 £ (— tt, tt], then 
either 9 = tt and, by (O, the sequence |f„| tends to oo, or G (— 7r,7r) and f„ converges to 
tg. But then by continuity, Z(tn) actually converges to Z(tg), which, as argued above, only 
coincides with e*^ when 9 — Q, that is, tg = 1/8. In this case, Z{tn) 1. 

■ 

Let us now go back to the evaluation of the tail distribution function of Nn via the 
integral Ij22|l . We choose a contour C — Cn that depends on n and consists of two parts ci^"* 
and Cn'' (see Figure EJ: 

• Cn^ is an arc of radius r„/8 = (1 + log^ n/n)/8, centered at the origin; note that its 
radius tends to 1/8 as n goes to infinity, 

(2) 

• C„ is a Hankel contour around 1/8, at distance l/(8n) of the real axis, which 
meets Cn^ at both ends; this contour shrinks around 1/8 as n goes to infinity; more 

(2) 

precisely, as t runs along C,\ , the variable z defined by 
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Figure 6. The integration contour 



runs over the truncated Hankel contour Hn shown on the right of Figured 

Un ^{x-i, xe [Q,Xn]} U {-e'\ e [-7r/2,7r/2]} U {x + i, xe [0,x„]} 

where (1 + Xn/nf' + l/n^ = r^, so that a;„ < log^ n and a;„ — log^ n + 0(l/n). 
We denote by z„ = a;„ + i the top right end of 7i„. This point tends to infinity as n does. 

The integral l|22|l on C = C„ is the sum of the contributions of the contours ci^^ and ci^^ 
We shall see that the dominant contribution is that of C„ , because of the vicinity of the 
singularity at t = 1/8. 

Let us first bound carefully Z{t) for t G C„. Let i„ € C„ be such that 

|Z(t„)Hmax|Z(i)|. 

By Lemma im |Z(t„)| tends to 1 as n grows. Moreover, every accumulation point a of the 
sequence tn satisfies \a\ < 1/8 and |^(a)| ~ 1. This forces a = 1/8, and we conclude that 
tn 1/8. Write t„ = (1 - it„)/8. Then u„ 0, but \u„\ > 1/n. By 

= 1 - 2ui/^ (1 + 0(1)) . 

Let us write, for short, -Un = 1 — Z{tn). Then Vn ^ but 

k„| =2|7.„|i/4(l + o(l)) >n-i/4 (26) 




Figure 7. The Hankel contour H and its truncated version H- 
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for n large enough. Moreover, 

1 TT 

|argi^„| = -|arg(w„)| +o(l) < - + o{\), 

so that 

cos(argw„) ^ + 

Finally, 

|Z(i„)p = |1 - i;„|2 = 1 - 2|i;„| cos(argt;„) + |i;„|2 < 1 - y2K>„| (1 + o(l)) , 

that is, 



|z(i„)l < 1 - ^Kl (1 + 0(1)) < 1 - 



The latter inequality follows from (l26|l . and holds for n large enough. Finally, for t G C„, 

l-|Z(t)|> VV4. (27) 

Let us now consider the integral on the contour Cn^ . By Lemma ITTl the quantity 

(1 + Zf Z^+i {1 + Z + Z^) (1 - zf 
1 + ^2 

is uniformly bounded on this contour by some constant c, independant of n and t. Moreover, 

|1 - > 1 - >l-\Z\> \n-^/^ 

by Ij27|l . The same bound holds for the term 1 — Z^^^. Therefore the modulus of the 
contribution of Cn^ in the integral (I22|l is bounded by 

4c8"ni/2 r-" ^ 0(8"ni/2-iog») ^ o(8"/n'") (28) 

for any m > 0. 

Let us now study the contribution of the contour C,\ . As t varies along C„ , the variable 
z defined by i = (1 + z/n)/'& varies along the contour Tin- As n goes to infinity, this contour 
converges to the contour TL shown on the left side of Figure Let z e H. Then z €Hn for 
n large enough, \z\ < \zn\ log^ n, and, as n goes to infinity, the following approximations 
hold with error terms independent of z: 

Z{t) = l-2{-zf'Si-^/* + o{n-^''^\ogn^ 
l-Z{t) = 2{-zf''^n-^/^ (l + Oin-^'^y/k^)) 

/ X (29) 

Z{ty = cxp(-2A(-z)i/'') f 1 + 0(n"i/"*logn)j (recall j = [\n^'^\) 

= 8"+ie-^ (1 + 0(log'* n/n)) . 

Observe that, for z (z H, the real part of (— z)^/'* is bounded from below by a positive 
constant a. Hence 

|exp(-2A(-z)i/4)| ^ cxp(-2A3?(-z)i/'') < cxp(-2Aa), 
so that exp(— 2A(— z)^/'*) does not approach 1. This allows us to write 

T^W - l-exp(-2A(-z)V4) (l + Oin-^/Ho,n)) . 
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12') 

Hence, uniformly in t € C„ , we have 

(1 + zf Z^+i (1 + Z + Z2) (1 - zf 



u,{t)t- 



(1 + Z2) (1 - (1 - 

smh (A(-z)V4) 



with 8t = 1 + z/n. Let us now integrate this over C 



(2). 



(2) 



6.8" 

6.8" 

372 



sinh2(A(-z)i/4) 



+ 0(1) 



•H 



sinh^(A(-z)i/4) 



We now put together our estimates of the integrals on Cn^ (Eq. ll28ll ^ and d)T' and obtain 

6.8"??-3/2 / c ^ze- 

Ih sinh2(A(-z)i/4)- 



2i7r 



-dz + 0(1) 



Using Ij2()ll and the estimation C„ ~ 4"n ^^^/y/ir, this gives 

3 /■ \/—ze~ 



> A) 



-dz. 



iy^ Jn sinh2(A(-z)i/4)' 

The next step in our proof of TheoremElis to set v = (— z)^/^ in the above integral. As z 
runs on H, the variable v runs on the contour J of Figure |H1 and the corresponding integral 
is easily seen to coincide with the integral on the contour T defined in the statement of the 
theorem. This gives the first expression of G(A). 




0.5 ,1 1.5 2 2.5 



Figure 8. The contours T (two half lines) and J. 

We now want to express G{X) as a real integral. We first observe that the integration 
contour F can be replaced by its translated version 

Fo = {-re-*''/^ r e (00, 0]} U {re"/^ r E [0, 00)}. 

This parametrization of Fq by r splits the integral into two real integrals, and one finds: 



12 

G(A) = 



TT Jo Vsinh^(Are'''/4) sinh^(Are-«'^/4) 
48 1 - cos(V2Ar) cosh(V2Ar) 5 



TT Jo (cosh(V2Ar) - cos(V2Ar))2 



dr. 
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The expected expression of G(A) follows, upon setting u = \f2\r. 

3.2. The limit law and its density 

We now want to prove that G'(A) is the tail distribution function of a random variable. 
Since it is the limit of non-increasing functions, it is non-increasing. Its integral expressions 
show that it is a continuous, and even a differentiable function of A on (0, -f-oo). In order to 
conclude, we still need to prove that Thm. 14.1] 

lim G(A) = and lim G(A) = 1. 

A— >oo A— +0 

In order to prove the first statement, we use the second expression of G(A) given in the 
theorem. We note that the function 

1 — cos u cosh u 



(coshu — cosu)2 

is well-defined, bounded and continuous on [0, +oo). Moreover, as u goes to infinity, 

1 — cos u cosh u 



so that the integral 



(coshu — cosmJ 

1 — cos u cosh 



0(e-"), 



/ 



C'du 



(coshu — cosu)^ 

is convergent. The term 1/A^ in the expression of G(A) then implies the convergence of G(A) 
to as A — > oo. 

In order to study the Hmit of G(A) as A ^ 0+, we consider instead the first expression of 
G(A). Since x^/sinh^(x) is analytic in the disk of radius tt, with expansion 1 — x^/3 + 0(a;^), 
there exists a constant c such that for \v\ < 7r/(2A), 

1 1 1 



smh\Xv) A2w2 3 



< cA 



(30) 



Let us write 



r sinh (Aw) 



-dv 



r Vsinh^(Aw) X^v^ 



v^e" dv - 



V 

A2 



V 

"3 



dv. 



Recall the Hankel expression of the reciprocal of the Gamma function, valid for any s G C: 



1 



T{s) 



Consequently, 



v^e^ 



dv — 



= 0, 



2r(o) 



'dz = — 
ill 



v^e" dv 



dv. 



(31) 



2r(-i/2) 



and we can rewrite 
12 



G(A) 



V"" JT sinh (Aw) 



-dv = 1 



12 



1 



1 



i ] v'^e"''dv. 



IT \s\uh^{\v) A2i;2 

Let us cut the above integral into two parts, \v\ < 7r/(2A) and \v\ > 7r/(2A). The first part 
is easily seen to tend to as A does, thanks to (I3()|l . For the second part, we observe that 
for A|w| > 7r/2, 

1 1 1 

smh\Xv) " A2^^3 

4 

is bounded (by a constant independent of A and v), that the integral of v^e"" on T is 
absolutely convergent, and that the contour {w € F : |u| > 7r/(2A)} "shrinks to oo" as A ^ 0. 
We finally conclude that G(A) tends to 1 as A ^ 0. 
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Consequently, there exists a random variable N having distribution function 1 — G(A), 
and Nn converges in law to N. Since G is differentiable, N has a density with respect to the 
Lebesgue measure on M+, which is /(A) = — G"(A). The two expressions of G given in the 
theorem provide the two expressions of /. 

3.3. The moments of N 

Let us first prove that for all fc > 0, the tail distribution function of N satisfies 

G{X) = o(A-'=) as A ^ oo. (32) 

This is easily seen to imply the existence of moments of N of all orders. In order to prove 
the above bound, we write 

^ ' iX^ e^^-" - 1 

This is obtained from the first expression of G(A) using an integration by parts. Now, for 
A > and u S T, 

|g2A. _ > |g2A.| _ ^ ^ g2AK(») _ 1 > e2A _ 

4 

Prom this, and from the term e" in the integral, we conclude that there exists a constant c 
such that 

G(A) < 

The bound Ij32ll follows. This bounds also guarantees that for A; > 1, 

E(A^'=) = fc / X^~^G{X)dX. (33) 

"'0 

The generic case: fc > 3. Recall the following integral representations of the Riemann 
zeta function: for 9fi(s) > 1, 

Cis) = — -— / aw = — / ^ aw = — / ^ ay. 

' r(s)io 4r(s + l)yo sinh2(«;/2) T{s + l)J„ smh\y)" 

The second expression follows from the first one after an integration by parts. 
Let us now combine H33II with the first expression of G(A): 

E(iV'=) ^ r x'^-'dX I "'f dv. (34) 

Assume for the moment that we can exchange the order of integration (this will be justified 
later). Exchange the integrals, and replace the variable A by y/v, where y is a new variable: 

E(Arfc) [ v'-'^e^'dv f -^^^dy. 

W"" JvR+ sinh (y) 

For fc > 3, the function y i-+ y'^^^ / smh^ {y) is meromorphic on C, with poles at ifcyr for fc g Z 
and fc ^ 0. From this, and from the strong decay of this function as ^{y) —* oo, it follows 
that the integral on y is actually independent of the choice of w S F. In particular, it is equal 
to its value at w = 1, which is 

°° y'-' ^ _ 4F(fc)C(fc-i) 

smh\y) ^ 2^ 

as recalled above. The integral on v is then evaluated in terms of the Gamma function 
using Ij31ll . and the expected expression of E(A^'^) follows. 

It remains to justify the exchange of integrals in (I34|l . Observe that 

I sinh(y)| = |e« - £"^'1/2 > (le^l - |e-^|) /2 = sinh(5R(y)), 
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SO that for v GT, 



1 



< 



1 



|sinh^(Aw)| ~ sinh2(A)' 

Moreover, the integral of v^e" along T is absolutely convergent, and so is the integral of 
A*^"^/ sinh^(A) over It follows that the integral (PIJ, once converted into two real 

integrals, is absolutely convergent, so that the integrals can be exchanged. 

The case fc = 1. We cannot apply exactly the same procedure as above, because the integral 
of 1/ sinh^(A) over is divergent. However, in view of ll3T|l . we can write 



G(A) 



12 



sinh^(Aw) A2i;2 



Also, replacing F by Tq in the latter integral does not change its value. The technique is 
then the same as above: 



E(iV) 



12 
12 



dX 



1 



1 



sinh^(Aw) A2t 



dv 



(35) 



sinh (y) y 



i ) 



(assuming we can change the order of integration). Again, the integral on y is independent 
of V, and equal to 

1 1 





'I 2 


oo 


)dy = 


y 


= -1 





,sinh^(y) y'^ , 
Using again (I31|l to evaluate the integral on v, one finds 



E{N) = 



r(-l/4) 2r(3/4)' 

In order to justify the exchange of integrals in (|35ll . we wish to prove that ll35)l is abso- 
lutely convergent. In order to do so, we split the integral over Fq into two real integrals, 
corresponding respectively to t; = re*^/^ and v = re~"/^. We are thus led to prove that 

1 1 



dX 



dr 



sinh^iXre'"/^) iX'^r'^ 

is absolutely convergent (and a similar result when i is replaced by —i). But we can exchange 
the order of integration in this integral of positive functions. Doing so, and setting X ~ y/r 
as above, proves that this integral is finite. 

The case k = 2. Let us start from another expression of G(A), obtained by writing v = w/X: 



G(A) 



12 



iy/nX^ 7;^P sinh^(w) 



12 



If 



v^A6 7p sinh^(w) 



,^'/''dw. 



The second expression follows from the analyticity properties of the integrand. Now, take 
e > 0, and let us evaluate 

AG(A)dA = / T^dA / — fl—e'^'/^^dw 



12 



sinh (ui) 



(36) 



w 



Vtt Jr sinh^(u' 



-dw 



A5 



-dX 



12 



Vtt Jp sinh^(w) 



w 



sinh (w) 



4 / 4 
j"' „ 



dw 



1 dw. 
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The absolute convergence of integrals that legitimates the exchange of integrals in Ij^ffiH is, 
this time, obvious (thanks to the fact that A > e). Now, the analyticity of the function 
w 1-^ w/sinh (w) for 5ft(w) > 0, and its strong decay as di{w) oo, imply that 



r sinh^(ii;) 



Hence 



i^TT 7r sinh^(it;) Jr sinh^(eu) 



XG{X)dX = ^ j ^ ^ — -dw = —f= I ^ 2, — 



3 f^_Ldv + oil) = ^ (by TO 



Now, observe that 

/oo 
XG{X)dX = E{N^lN>e) - e'^G{e). 

The announced expression of the second moment of N follows. 

3.4. Convergence of the moments of iV„ 

In this section, we prove that the moments of Nn — Mnjv}!^ converge to the correspond- 
ing moments of N . In order to do so, we first express E(M^) as the coefficient of i" in a 
certain series. Then, we apply the general consequences of the analysis of singularities: if 
this series is regular enough (with a precise meaning of regular) , one can derive the asymp- 
totic behaviour of its coefficients from the singular behaviour of the series near its dominant 
singularities |17) . 

Recall that the series JT,, given by Ij21|l . counts the trees that contain at least one label 
larger than j. Hence Uj-\ — Uj counts the trees having maximal label j. Also, note that 

Uj = V[Z^) - V{Z^+^), (37) 

where 

^ xZ(l + Z)(l-Z3) 

^ ' (l + Z2)(l-xZ2) ■ 

Consequently, for fc > 1, 

m-i!:) = ^ E^''n(^.-i - u,) = -^n {u + 1)^ - j') u,. m 

For k = 1, this gives 

2"C„E(Af„) = n ^ {V{Z^) ~ V{Z^+^)) = r] [V{1) + V{Z)) = H ^ ^ , '-. 

By Lemma im the latter series is analytic in C \ [1/8, oo). The generic consequences of the 
analysis of singularities apply: one can derive the asymptotic behaviour of the coefficients 
from the singular behaviour of the series Given that, when t 1/8, 

Z{1 + 2Z + 2Z'^) _ 5 



^^^2 = g - 6(1 - 8i)i/4 + 0(^/^^), 

the behaviour of the nth coefficient of this series is 

, „,Z(1 + 2Z + 2Z^) ^S^n-^/"^,, ,,,, 3 8"n-5/4^^ 

It remains to divide by 2"C„ ~ 8"n~^/^ /y/n to conclude that 
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which is also the first moment of N. 

Now, by combining the expression (13 711 of Uj and (I38|l . one obtains, for k >2, 

2"C„E(A/,t) - n I Vil) + (2^- - 1)V{Z) + ^ {U + If - f - (j - If + (j - 2)'^) F(Z^) 

(39) 

Observe that {j + if — j'^ — (j — 1)'^ + {j — 2f is a polynomial in j of degree k — 2 and 
leading coefficient 2k{k — 1). Let 

We are going to prove that, for ^ G N, 

f \^ 

I 2^r(€/4) 

Assume for the moment this is proved, and let us conclude about the limiting moments of 
iV„ = M„n-i/4. First, we observe that for j > 0, V{Z^) has (only) a fourth root singularity, 
so that the coefficient of in V{Z^) grows like 8"n~^/"', up to a multiplicative constant. 
This observation, combined with Ij39|l and the above asymptotics of a„(€), implies that the 
dominant term in the asymptotic behaviour of 2"C„E(M^) is that of 2k{k — l)a„(A; — 2). 
After normalizing by 2"C„n'^/^, this gives 

( 3V^ if A: = 2, 

E(M>-^/4)^J 24V¥fc!C(fc - 1) 

I 2'=r((fc 2)/4) ^f"-^- 

These limiting moments are exactly those of N . 

It remains to study the asymptotic behaviour of the numbers a„(£) (for i fixed, and n 
going to infinity). We have: 

^^^^^ - Z(l + Z2) j^^'T^ 
(1 + Z)(1-Z3) ''' .„ 
Z(l + Z2) ^ 

^ ' j>l,m>l 

(1 + Z)(1-Z3) 



Z{l + Z^) ^ 

^ ' N>1 



where 



The function 



3\N 



N>1 

is easily seen to have radius of convergence 1. Moreover, as z tends to 1 in such a way 
|arg(l-z)| <0< V2, 



Diiz) 



^ log ( -— ) if £ = 0, 



1 - Z V 1 - 2 

(l-z)^+i 



if ^> 1 
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(this can be obtained using a Mellin transform \W\ 115)1. The above expression of Ai{t), 
combined with Lemma ITTl and these properties of Di{z), shows that Ai{t) is analytic in the 
domain 2? = C \ [1/8, oo). Moreover, since | arg(l — Z)\ < 7r/4 + o(l) a,s t 1/8 in V, we 
can use the above estimates of Di{z). This gives 



Aiit) 



-31og2 + ^log(^^ 
3^!C(€+1) 



8< 



if - 



0, 



if ^ > 1. 



2^(1 - sty/^ 

The generic results derived from the analysis of singularities apply, and give the asymptotic 
behaviour ijlfUl of the numbers a„(£). This concludes the proof of Theorem El 

3.5. The supremum of the support of the ise 

Let us finally prove PropositionEH The following argument requires a detour via discrete 
snakes and Brownian snakes. We refer to [2ll [231 1201 for definitions and notation^. In 
particular, we use the following integral representation of the random measure ^jse^ for any 
continuous bounded function g on M, 



g{y)dnisEiy) 



g{r{t))dt 



(41) 



where r(.) is a random process, continuous on [0, 1], called the head of the Brownian snake. 
In other words, /iise is the occupation measure of the process r. (Again, the definition of r 
varies from one paper to the other. The above formula fixes our normalization of r.) 

The random variable Nn = Mnn~'^/^ coincides with max(r„), where r„ is the (normal- 
ized) head of the discrete snake associated with our tree family. The random process \/2r„ 
converges weakly to r, the head of the Brownian snake (23]. Since max is a continuous func- 
tional on C[0, 1], this implies that = '\/2max(r„) converges in distribution to max(r). 
Thus max(r) has density f{\/y/2)/^/2, where / is defined in TheoremEl 



A A + e 



A- 




A - e A A + e 



Figure 9. The functions g\,e and /iA,e. 

It remains to prove that max(r) is equal (in distribution) to Nig-^, the supremum of the 
support of the ISE. Let A G M and e > 0. Let be the function plotted on the left-hand 
side of Figure [HI We have 

iVisE < A ^ise(-oo, A] = 1 -^=+ / fx,i{y)dfiisEiy) = 1. 

Jm 

Thanks to (QlJ, this gives 

iVisE<A^ / h.Ar{t))dt = 1 



We warn the reader that normalizations change from one paper to another. 
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Taking probabilities yields to 

P(^iSE < A) = P / fx,Mt))dt = 1 ) = P(max(r) < A), 
since r is almost surely continuous. 







4. A LOCAL LIMIT LAW 

For j G Z, let Xn{j) denote the random variable equal to the number of nodes having 
label j in a random tree of 7o „. This quantity is related to the series Sj{t,u) studied in 
Proposition 01 In particular, 



Also, observe that 

^«(j)-O^A/„ <j, 

where M„ is the largest label, studied in the previous section. Let us define a normalized 
version of X„(j) by 

Let A e M. The aim of this section is to prove that F„([An^/'*J) converges in distribution, 
as n goes to infinity, to a random variable Y{X) that we describe by its Laplace transform. 
This is achieved in Theorem^lbelow, but we first want to present two consequences of this 
theorem, which have a simpler formulation. The first consequence deals with the case A = 0. 
Recall that, up to a normalization by n'^/'*, the random variable Yn{0) gives the number of 
nodes labelled in a tree rooted at 0. 

Proposition 12 (The number of nodes labelled 0). As n goes to infinity, the random 
variable 3y„(0)/^/2 converges in distribution to T where T follows a unilateral stable 
law of parameter 2/3. The convergence of the moments holds as well: for k >0, 



E(y„(o) 




^ r(i + 3fc/4) 
r(i + /c/2) 




This proposition will be proved in Section IT^ I am indebted to Alain Rouault, who rec- 
ognized that the above moments were related to T. Recall that T is given by its Laplace 
transform: 

Eie-"^) = e""'''" for a > 0. 

The second consequence of Theorem ^1 is an explicit expansion in A of the limiting first 
moment of Yn{j). 

Proposition 13 (The first moment). Let A £ K. Denote j = [An^/^J . Then, as n goes 
to infinity, 

Star CO. <=l±i)^r(^ 

^ m>0 ■ ^ 

This function of \ is plotted on Figure nUl 

Similar, but more and more complicated expressions may be written for the next moments of 
Y{X). This proposition will be proved in Section lT^ Let us, finally, state our main theorem, 
from which the two above propositions derive. 
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Figure 10. The average number of nodes labelled [A-y/nJ in a tree of size 
n, when n oo. 

Theorem 14 (A local limit law). Let A > 0. The sequence Yn{[Xn^^'^\) converges in 
distribution to a non-negative random variable Y{X) whose Laplace transform is given, for 



More precisely, the Laplace transform o/ y„( [An^/**] ) converges pointwise to L{X,-) on the 
interval (— 4/a/3, 4/a/3). The convergence of moments holds as well. 

It is believed (or known?) that the random measure ^ise is almost surely absolutely 
continuous with respect to the Lebesgue measure on M. Eq. Q leads us to the following 
conjecture. 

Conjecture 15 (The density of the ISE). There exists a random continuous process 
/ise(A), defined for A e M, such that //ise = fisELeb, where Leb denotes the Lebesgue measure 
on M. Moreover, /ise(A) satisfies 



where the law ofY{X) is given in Theorem,\14\ 
Comments 

1. The Hmit random variable Y{X) equals with a positive probability as soon as A > 0. 
Indeed, by the portmanteau Theorem |14[ Thm. 11.1.1], 



\a\ < 4/V3, by 




where 




(42) 




P(r(A) = 0) > limsupP(y„([An^/4j) = 0) = limsupP(A/„ < [Xn^ 



/4J). 



But, by Theorem El 



P(M„ < LAn^/^'j) ^ 1 - G(A) > 0. 
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2. Let US add a few words on the series A{x) defined by Ij42|l . in order to convince ourselves 
that the integral giving L(A, a) is well-defined. Clearly, the expansion of A{x) at x = 
has non-negative coefficients. Looking at the discriminant of the equation that defines A 
shows that A has radius of convergence at least 4/a/3. Moreover, it is easy to prove that 
A(4/V3) = 2 - V3 = 0.26 . . . Consequently, \A{x)\ is bounded by 2 - for \x\ < 4/V3. 
Since \v\ > 1 for w e F, the modulus of A{a/v^) is bounded from above by 2— a/S. Moreover, 
> 1, so that [e-^^"] < e'^^ < 1. Hence 

(l+A(a/w3)e-2^'')2 

is uniformly bounded on F, and L{X, a) is well-defined. 

Note that the series A{x) admits the following closed form expression: 

A(x) = p 1. (43) 

l + _2^cos(2i2£££(_£^) 

This can be checked by proving that this expression satisfies Ij42|l and the initial condition 
A{0) = 0. 

4.1. Proof of Theorem 1141 

Let A > and j ~ [Ani/4J . Let us first express the Laplace transform of Yn{j) in terms 
of the generating functions Sj{t, u) of Proposition ^ 

,aY„U)\ = E fe-"-'^'^^U)) = (44) 
/V / 2^C^ 

Again, we will evaluate this Laplace transform thanks to the analysis of singularities |17j . 
We wish to use again the integration contour C„ of Figure This requires to prove that 
Sj{t,u) is analytic in a neigborhood of this contour (for n large and u = e"" ' ). This is 
guaranteed by the following lemma. This lemma naturally includes some properties of the 
series ^ involved in the product form (jllll of Sj . We denote by Jn the part of the complex 
plane enclosed by C„ (including C„ itself) . 

Lemma 16 (Analytic properties of /i and Sj). Let a be a real number such that \a\ < 
4/-\/3. Then there exists e > such that for n large enough, the series /i(t,M„), with u„ = 
, is analytic in the domain 

^ {t:\t- 1/8| > l/((8 + e)n)} \ [1/8, +(»). 

In particular, ^(i,u„) is analytic in a neighborhood o/X„. Its modulus in X„ is smaller 
than a, for some a < I independent of a and n. The series Sj{t,Un) is also analytic in a 
neighborhood o/X„. 

Proof. The lemma is clear if a = 0: in this case, u„ = 1, the series /i(t, u„) vanishes, and 
the series Sj reduces to the size generating function of labelled trees, namely T, which is 
analytic in C \ [1/8, oo). We now assume that a ^ and |a| < 4/V3. This guarantees that 
A{a) is well-defined, where the series A is defined in Theorem ITU 

Let us first study the singularities of the series p = fi[z, u) defined as the unique formal 
power series in z satisfying 

_^ {l + Z^){l+rLz){\+-^Z^){\+JjLZ^) 

^(1 + Z)(1 + Z + Z2)(1 -2)3(1 _/2z2)- 

Note that has polynomial coefficients in u, and vanishes when u = I. Assume that it is a 
fixed real number close to, but different from, 1. Recall that, as all algebraic formal power 
series, p,{t, u) has a positive radius of convergence. Let us perform a classical analysis to 
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detect its possible singularities. These singularities are found in the union of two sets Si and 
S2: 

• Si is the set of non-zero roots of the dominant coefficient of the equation defining p,. 
That is, Si = {±i} , 

• ^2 is the set of the roots of the discriminant of the equation defining p,. For u ^ 1 + x 
and X small, these roots are found to be 

z = ±l, z = -l + 0{x), z = e±2''^/3 _^ 0(x), z^l+ujl2^/^x^/^ + Oi\x{^/^), 

where cu satisfies ui^ = 1. (The term uj allows us to write loosely x-'-/^ without saying which 
determination of the cubic root we take.) 

Observe that the moduli of all these "candidates for singularities" go to 1 as x goes to 0. 
Now the series ^ = fJ.(t,u) involved in the expression Ijllll of Sj{t,u) satisfies 

fj,{t,u) = fl{Z{t),u) 

where Z{t) is defined by In other words, we could have defined the series ft by 



H{z,u) = 11 



(1 



Recall that Z is analytic in the domain P = C \ [1/8, 00). Take u = Un ~ e°" ' = 1 + x, 
with X = an~^/^[\ + o(l)). By Lemma ITTl for n large, the only values of Si U ^2 that may 
be reached by Z{t), for t G T), are of the form 

In view of 101, these values of Z{t) are reached for 

8 128 n ^ ' 

Since \a\ < 4/-\/3, these values of t are at distance less than l/((8 + e)n) of 1/8, for some 
e > 0, and hence outside the domain £„. Consequently, u„) is analytic inside £„. 

We now want to bound ^(t, w„) inside X„. Let t„ G X„ be such that 

lM(tn,Wn)| = max|^(i, u„)|. 

In particular, > |/i((l — l/n)/8, Un)\. In order to evaluate the latter quantity, note 

that Z((l — l/7i)/8) = 1 — 2n^^/'' + 0(n^^/^). Thanks to the closed form expression of ^ given 
in PropositionEl and to the expression of the series A, we see that ^/n)/8,Un) —* 

A{a). Since a 7^ 0, A{a) ^ 0, and for n large enough, 

Htn,u„)\ > |^((l-l/n)/8,M„)| - \A{a)\ + o{l) > 0. (45) 

Recall that all the sets X„ are included in a ball of finite radius centered at the origin. Let 
a be an accumulation point of the sequence t„. Then \a\ < 1/8. 

Assume first that a ^ 1/8. Then there exists N such that a is in for all n > N, that 
is, in the analyticity domain of ij{-,Un). Let tn^,tn^, . . . converge to a. By continuity of /x 
in t and u, we have 

M(ini,Mn,) /i(ck, 1) = 0. 

This contradicts Hence the only accumulation point of t„ is 1/8, and t„ converges to 

1/8. Let us thus write 

1 / X,-, 

t„ = ^ 1-- 
8 v n 

We have Xn = o{n), but also \xn\ > 1 since t„ belongs to I„. We wish to estimate /i(i„, m„). 
From the singular behaviour of Z (Lemma II l|l . we derive 

Z(t„) = l-2(^)'^Vo('(^''^' 
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Moreover, 
This gives 



1/4N 

_ I 1 4- n I I ^ ' 

If the sequence Xn was unbounded, then there would exist a subsequence x„. converging to 
infinity. Then — 1)/(1 ~ Z)^ would tend to 0. The closed form expression of /i given 
in Proposition implies that ^(tm^Un.) would tend to 0, contradicting Ij45|l . Hence the 
sequence Xn is bounded, and one derives from the explicit expressions of fi and A that 

M(in,u„) = A{ax-^/'^) +o{l). 

Since A is bounded by 2 — V3 inside its disk of convergence, |/i(t„, m„)| is certainly smaller 
than some a for a < 1 and n large enough. This concludes the proof of the second statement 
of Lemma nm 

By continuity of ^{t,Un), this function of t is still bounded by 1 (in modulus) is a neigh- 
borhood of X„. Recall also that the modulus of Z(t) never reaches 1 for t G C \ [1/8, oo). 
The form lfTT|l then implies that Sj{t,Un) is an analytic function of t in a neigbourhood of 

T 

■ 

Let us now go back to the expression ijliji of the Laplace transform of Yn{j). Thanks to 
the lemma we have just proved, we can use the Cauchy formula to extract the coefficient of 
in Sj{t,Un). We use the following expression of Sj: 

' (l + /iZJ+l)(l + /iZJ+3) ' 

which is easily derived from (lllll . Thus 



Again, we split the contour C„ into two parts CrP and Cn \ shown in Figure El As in the 
proof of Theorem El the contribution of Cn^ is easily seen to be o(8"/n™) for all to > 0, 
thanks to the results of Lemmas II II and [T?)l On Cn\ one has 



8 V 

where z lies in the truncated Hankel contour Conversely, let z G H.. Then z G Tin 

for n large enough, and, in addition to the estimations ll29|l already used in the proof of 
Theorem ini one finds 

^l{t, w„) = ^(a(-z)-3/4)(i + o(l)), (46) 
where A{x) is the series defined by lll2)l . After a few reductions, one finally obtains 



12.8"n-3/2 r A(a(-z)-3/4)cxp(-2A(-z)i/4) 



[t \b-j{t,Un) ^2 6„H : / - , , .^/.\ f n\i \Mi\\7 d'Z+o(8 n ' ). 



(1 + A(a(-z)-3/4) cxp(-2A(-z)i/4))2 

It remains to normalize by 2"C„ = W^n^'^^'^ j ^ , and then to set v — {—zy/'^ to obtain the 
expected expression for the Hmit of the Laplace transform of Yn{j), with j = [An^/'*J . 

The limit Laplace transform L{X, a) is clearly continuous at a = 0, and equals 1 at this 
point. A version of Levy's continuity theorem ^1 Thm. 9.8.2] adapted to Laplace transforms 
implies that the sequence Yn{j) converges in distribution to a limit random variable Y{X) 
having Laplace transform L(A, •). 
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From the convergence of the Laplace transform in a neighbourhood of the origin, it is 
easy to derive that for every k, the sequence of random variables Yn{j)'' is uniformly inte- 
grable. But then the convergence in distribution implies the convergence of the moments (4) 
Thm. 5.4]. This concludes the proof of Theorem 1141 



4.2. Proof of Proposition fT2l 

When A = 0, the limiting Laplace transform reduces to 

L(0, a) = 1 + —= / TTT'T^'"^ dv = l + \ / sn ''' e ^« 



(l + A(a/f3))2 ' iV^7p l + x(a/t'3) 

where x(a^) is the unique series in x satisfying 

x = ^(i + x)'/^. 

The Lagrange inversion formula |25[ p. 38] gives, for fc > 1, 



„fei 



X{x) _ 1 r(3fc/2 ~ 1) 



' 'l + xW 6*^- fc!r(fc/2) 

Consequently, 

r/n N 1 12 /"^ 1 r(3fc/2-l) 5_3fe 

The convergence is absolute, so that we can exchange the sum and the integral: 

r/n N . 12 ^ 1 r(3fc/2-l) /■ 5_3fc ,„4^ 

Using (EU, and picking the coefficient of , we find that the /cth moment of the random 
variable F(0) is 

, , r(3fc/2 - 1) 

^ w ; gfc_i r(/c/2)r((3fc-2)/4)' 

The duplication formula, 

22«-ir(s)r(s + i/2) = 0Fr(2s), 

appHed to s = (3fc — 2)/4, finally gives 

Since F(0) has a Laplace transform, it is uniquely determined by its moments Thm. 30.1]. 
But 

_ r(i + 3fc/4) 
r(i + fc/2) 

is known to be the fcth moment of T^^/-^, where T follows a unilateral stable law of parameter 
2/3 (see p. ill]). Proposition [H follows. 
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4.3. Proof of Proposition fisl 

We have derived above the moments of Y{0) from the expression of its Laplace transform. 
This extends to the moments of Y{X), for A > 0: for A: > 1, 

Since A(x) = x/2A + 0{x^), the case /s = 1 of the above identity reads 

E(y(A)) = / e-^^'^v^e^'dv. 



iV^^ Jv 

In the above expression, expand the exponential as a series. The convergence of the sum 
and integral is absolute, so that one can exchange them. This gives: 



E,y(A), . 4 j: / 

T7l>0 



I' e 

r 



Using Ij31ll (which is valid for any s with the convention \/T{—n) = for n G N), this can 
be rewritten as 

E r A) =^fn} — TFTT^ — , = ^ > ^ — - — cos ^ 

^ " ^mr((l~m)/4) x/tt^to V4 / V 4 

m>0 / V \ / \ 

The last equality follows from the complement formula, 

5. A global limit law, and the distribution function of the ISE 

In Section 21 we have derived from Proposition 21 some local limit results; for instance, 
a limit law for X„(0)/n'^/*, the (normalized) number of nodes labelled 0. In this section, 
we proceed with a similar study, which aims at deriving from Proposition a global limit 
result — in particular, the limit law of X+(0)/n, the normalized number of nodes having a 
non-negative label. The technique is copied on Sectional and we do not give all the details. 

For j G Z, let X^{j) denote the random variable equal to the number of nodes having 
label at least j in a random tree of 7o_„. Let us define a normalized version of X^(j) by 

r-(,) = SO). 



These quantities are related to the series Rj{t,u) studied in Proposition In particular. 

We extend the definition of X+ and Y^f to real values in a natural way by setting X+ (a;) = 
X+([a;]) and Y^(x) = y,+ ([a;]). Let A > 0. The aim of this section is to prove that 
Y^i^X-n}/^) converges in distribution, as n goes to infinity, to a random variable ^'''(A) that 
we describe by its Laplace transform. This is achieved in Theorem E| below, but we first 
want to present two consequences of this theorem, which have a simpler formulation. The 
first consequence is a striking limit law for y^(0). Recall that, up to a normalization by rt, 
this random variable gives the number of nodes having a non-negative label in a tree rooted 
at 0. 

Proposition 17 (The number of non-negative nodes). As n goes to infinity, the ran- 
dom variable Y^{0) converges in law to the uniform distribution on [0,1]. 

This proposition will be proved in Section EH The second consequence of Theorem E| is an 
explicit expansion in A of the limiting first moment of Y^{\n^/^). 
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Proposition 18 (The first moment). Let A > 0. Then, as n goes to infinity, 

■=(-"^'-"'))-3^s:^»(x)^ 

^ m>0 

This proposition will be proved in Section lOI 

Let us, finally, state our main theorem, from which the two above propositions derive. 

Theorem 19 (A global limit law). Let A > O.The sequence {Xn^^^) converges in 
distribution to a random variable Y^{X) whose Laplace transform is given, for \a\ < 1, by 



where 




G(A,a) = 1 H ^ / ^ ; , — -^--v'e" dv, 

and the integral is taken over 

r = {1 - ^e-"/^ t e (oo, o]} u {i + ^e-*'^/^ t e [o, oo)}. 

Moreover, the Laplace transform ofY^{Xn^/^) converges pointwise to G{X, •) on the interval 
(— 1, 1). The convergence of moments holds as well. 

The sequence 1^(— An^/*) converges in distribution to the random variable 1 — y+(A). 

This theorem will be proved in the next subsection. In view of the following proposition, 
it tells us about the law of the distribution function of the ISE. 

Proposition 20 (The tail distribution function of the ISE). Let 5ise(A) = 
A*ise(A, +oo) denote the tail distribution function of the ISE. Then for A > 0, 

.9ise(A) = Y+{X/V2), 

where the law of the variable ^^(A) is given in Theorem,\19[ In particular, 5ise(0) is uni- 
formly distributed on [0,1]. The random variable 5ise(— A) has the same distribution as 
1-Y+{X/V2). 

Comments 

1. The law of 5ise(0) was already given by Aldous Eq. (12)]. 

2. Let us add a few words on the series B and D to convince ourselves that the integral 
giving G(A, a) is well-defined as long as \a\ < 1. Let E{x) = 1 — D{x). Then E admits the 
following expansion: 




n>l 

where C„ is the nth Catalan number. Similarly, 

1 - 2 ^ G„_i4-"a;", 

n>l 

and these two identities imply that E{x) has non-negative coefficients. Moreover, its radius 
of convergence is easily seen to be 1, so that |i?(a;)| < E{1) = 1 — l/\/2for |a:;| < 1. Moreover, 
expressing B in terms of E gives: 

E{1 - 2E) 



B 



{2- E){i~2Ey 
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which shows that B{x) is also analytic for |a;| < 1 and satisfies in this domain 

_B(1)(1 + 2^(1)^ /- 

For V er,\v\>l and ^{v) > 1. This implies that 

B(a/'y4)e-2^'' 



(1 + B(a/w4)e-2>'')2 
is uniformly bounded on F, and G(A, a) is well-defined. 

5.1. Proof of Theorem 1191 

Let j = [An^/"*] . Given that the product forms for the series Sj and Rj are very similar, 
it is not surprising that we use an approach copied on that of the previous section. We start 
from 

with Un = e°/". For technical reasons, we choose to modify slightly the integration contour 
of Figure The Hankel part of this contour, which was lying at distance 1/8 of the real 
axis, is now moved a bit further, at distance 1/6 of the real axis. More precisely, the new 

contour C„ consists of two parts C^^^ and C^^ such that 

• C„ is an arc of radius (1 + log n/n)/8, centered at the origin; 

(2) 

• C„ is a Hankel contour around 1/8, at distance l/(6n) of the real axis, which meets 

at both ends. 

We first need to prove that the series Rj(t,Un) is analytic in a neighborhood of X„, the 
region lying inside the integration contour C„. The following lemma is the counterpart of 
Lemma llfil 

Lemma 21 (Analytic properties of v and Rj). Let a be a real number such that \a\ < 1. 
Then i/{t,Un) is analytic in a neighborhood o/X„. Its modulus in X„ is smaller than a, for 
some a < 1 independent of a and n. The series Rj{t,Un) is also analytic in a neighborhood 

of In- 

Proof. Again, the lemma is obvious if a = 0. We thus assume a ^ and \a\ < 1. 
Let us first study the singularities of the series P = v{z,u) defined by 

According to Proposition Q v is a formal power series in z with polynomial coefficients in u, 
and by 10, one has: 

v{t, u) — v{Z{t), u). 

In the course of the proof of Proposition we have obtained a polynomial equation 
P{u,Z,u) = 0, of degree 4 in relating v{t,u), Z{t) and the variable u. This equation 
is not written in the paper (it is a bit too big), but it can be easily obtained using the 
expression of v given in Proposition |H1 By definition of P, we have P{v, z, u) = 0. 

Assume that m is a fixed real number close to 1. That is, u = 1 + cc, with x small. In 
order to study the singularities of we look again at the zeroes of the leading coefficient of 
P, and at the zeroes of its discriminant. This gives several candidates for the singularities of 
D{z, u), which we classify in three series according to their behaviour when x is small. First, 
some candidates tend to a limit that is different from 1, 

z = -l, z^±i, z = e±2W3^ z = e±2W3 ^ 0(2,)^ z = -l + 0{x). 
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Then, some candidates tend to 1 and lie at distance at most \x\^^^ of 1 (up to a multiplicative 
constant): 

where w is a fourth root of unity and c is in the set {0,16,64/3,-16/3}. Finally, some 
candidates tend to 1 but lie further away from 1 (more precisely, at distance |a::|^/^): 

z = l + 2e^"/Vxi/^ + 0(|a;|i/3), 

where uj' is a sixth root of unity. 

Let us now consider i/(t, u) = D{Z{t),u) with u — Un ~ e"/" = 1 + x, where x = 
a/n(l +o(l)). Recall that Z is analytic in C\ [1/8, oo). By Lemma ITTl the series Z{t) never 
approaches any root of unity different from 1. Hence for n large enough, Z[t) never reaches 
any of the candidates z of the first series. 

The candidates of the second series are of the form 

z = l+uj{ac/nf'^ + 0{n-~^'^) 

for some constant c, with |c| < 64/3, depending on the candidate. By (O, Z{t) may only 
reach these values for 

Since \a\ < 1, there exists e > such that these values lie at distance less that l/((6 + e)n) 
of 1/8, that is, outside a neighborhood of the domain X„. 

The candidates of the third series are more worrying: Z{t) may reach them for 

t^l-'^ia/nf' + Oin-^'^), (49) 

where w" is a cubic root of unity, and these values may lie inside X„. If a > and uj" = 
g±2t7r/3^ or if a < and w" = e^"^/^, the modulus of the above value of t is found to be 
1/8(1 + crr"^/^ + o(ri^^/^)), for some positive constant c: this is larger than the radius of the 
contour C„, which implies that t lies outside a neighborhood of I„. However, if a > and 
w" = 1, or if a < and lo" = 1 or e"^*'^/'^, the above value of t lies definitely inside X„. Its 
modulus is 1/8(1 — crr'^l'^ + o(n~^/'^)), for some positive constant c. 

In order to rule out the possibility that iy{t,Un) has such a singularity, we are going to 
prove, by having a close look at the expression of ly given in Proposition that the radius of 
convergence of z^(t,u„) is at least 1/8 — 0(l/n). We use below the notation of Proposition |H1 

Clearly, the series V{t,Un) has radius of convergence min(l/8, l/(8u„)). In particular, 
this radius is at least pn ■— 1/(8(1 + \x\)) (with m„ = 1 + x). Moreover, the series V admits 
the following expansion 



1 / / 8tx \ I ^ „ f 2tx 



/ n>l 



4 \ V 1 - 8t / 2 ^ V 1 - 8i 

\ / n>l ^ 

where C„ is the nth Catalan number. This shows that V{t,l + |a;|) is a series in t with 
positive coefficients and that for all t such that \t\ < p„, 

\V{t, l + x)\< V{\t\,l + \x\)<V (^(YTnT' ^ ^ " ^' 
The next step is to prove that A(t, u„) never vanishes for \t\ < pn- Indeed, 

A = (1 - F)2 - 4W^F^ 

where W = W{t) is the formal power series in t defined by (I24|l . This series has radius 1/8, 
and non-negative coefficients. Hence for all t such that \t\ < 1/8, one has |W(t)| < W{l/8) — 
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p2 



1/4. Consequently, for \t\ < p„, 

m, 1 + x)i > (1 - \vit, 1 + x)\r - mmvit, i + .t)^ > (^i - ^ j -j^^i 

Hence A(t,M„) does not vanish in the centered disk of radius p„. It follows that the series 
P{t,Un) is analytic inside this disk. 

According to the expression of v given in Proposition|Hl the series u(t, Un) is meromorphic 
for |t| < pn- The final question we need to answer is whether v has poles in this disk, and 
where. Returning to the polynomial P such that P{v, Z,u) =0 shows that this can only 
happen if the coefficient of i^'^ in this polynomial vanishes. But this can only occur if z = Z{t) 
has one of the following forms: 

z = ±i, z = e±2W3^_o(2:), z = -l + 0{x), z = 1 + uj{64x/3)^/^ + 0{x^^^). 

As argued above, only the last value of z is likely to be reached by Z{t), and this may only 
occur if 

t=l-lia/n) + Oin-'^'). 

Consequently, the radius of v(t,Un) is at least 1/8 — 0{\/n), and this proves that the val- 
ues Ij49ll that have been shown to lie in the centered disk of radius 1/8, are not, after all, 
singularities of v{t, Un). This completes our proof that i'{t, u„) is analytic in a neighborhood 

of In. 

We now want to bound iy{t, u„) inside X„. From now on, we can walk safely in the steps 
of the proof of Lemma [161 Let t„ € X„ be such that 

Witn,Un)\ = max u„)|. 

t6l„ 

We first give a lower bound for this quantity, by estimating v{t,Un) for t = 1/8 — l/(6n). 
This is easily done by combining the closed form expressions of v (Proposition |Hl and B 
(Theorem EJ. One obtains: 

W{tn, un)\ > \p{l/8 - l/(6n), u„)| = \B{3a/4)\ + o(l) > 0. 

This lower bound is then used to rule out the possibility that the sequence t„ has an accu- 
mulation point different from 1/8. Thus t„ converges to 1/8, and one can write 

1 / Xr 

t„ = - 1 - — 



n 

We have Xn = o(n), but also |x„| > 4/3 since i„ belongs to X„. We want to estimate 
u{tn,Un)- Since 

"r \ 1/4 / /T N 1/2^ 

and 

one has 



u„ - 1 = a/n (1 + 0{l/n)) , 
u- 1 



(1 - Z)4 16a 

The closed form expressions of u and B imply that the sequence a;„ is bounded and 

v{tn,Un) = B{a/Xn) + o(l). 

Since B is bounded by 0.12 inside its disk of convergence, \i^{tn, u„)\ is certainly smaller than 
some a for a < 1 and n large enough. This concludes the proof of the second statement of 
Lemma 121 

By continuity of j/(<, ii„), this function of t is still bounded by 1 (in modulus) is a neigh- 
borhood of T„. Recall also that the modulus of Z{t) never reaches 1 for < e C \ [1/8, cx)). 
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The form Ijl6|l then impHes that Rj(t,Un) is an analytic function of t in a neigborhood of 
T 



Once this rather painful lemma is at last established, the rest of the proof of Theorem 1 191 
copies the end of the proof of Theorem with Sj , ^ and A respectively replaced by Rj , v 
and B. The counterpart of Ij46|l is 

v{t,Un)^B{-a/z){l + o{l)). 

Recall that the Hankel part of the contour C„ is now at distance l/(6n) of the real axis. 
Hence, when n goes to infinity, one finds 

74/3^ (l + S(-a/z)cxp(-2A(-z)i/4))2 
After normalizing by 2"C„ and setting v = (— z)^/^, this gives 

but the analyticity properties of the integrand allow us to replace the integration contour by 

r. 

5.2. Proof of Proposition [T3 

When A = 0, the limiting Laplace transform reduces to 

G(0, a) ^1 + —= / /-,,□./ 4NN2 f e dv = \ + -—= / — — , , v e dv 



where x(a;) is the unique formal power series in x satisfying 

The Lagrange inversion formula gives, for fc > 1, 

r^feiX(^)(3-x(a:)) 6 (2fc - 2)! 



' ' 1 + X(2;) 4*^ (fc- l)!(fc + 1)!' 

Consequently, 

1 (2fc-2)! 



^ -'^ k>\ 



4*^ (fc - l)!(fc+ 1)! 



The convergence is absolute, so that we can exchange the sum and the integral: 

../¥^4fc(fc-l)!(fc+l)! ir 



fc>i 



Using (EU, and picking the coefficient of a'^, we find that the fcth moment of the random 
variable i^^(O) is 



i0F4fc (fc- l)!(fc + 1) 2r(fc - 1/2) fc+r 

The unique distribution having its fcth moment equal to l/(fc + 1) is the uniform distribution 
on [0,1]. Proposition II 71 follows . 
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5.3. Proof of Proposition fisl 

We have derived above the moments of y+(0) from the expression of its Laplace transform. 
This extends to the moments of Y^{X), for A > 0: for fc > 1, 

Since B{x) ~ .t/48 + 0{x^), the case fc = 1 of the above identity gives 

E(r+(A)) = / e-^^^we^'dv. 



l\/TT 



r 



In the above expression, expand the exponential as a series. The convergence of the sum 
and integral is absolute, so that one can exchange them. This gives: 



E(r+(A)) = / v"^+'e^"dv 



Using l(3T|l . this can be rewritten as 



2 ^ m!r((2-m)/4) 2J^ ^ to! \ A V 4 

m>0 " " ' ^ m>0 ^ ^ 

The last equality follows from the complement formula (I47|l . 
5.4. The distribution function of the ISE 

Let us finally prove Proposition 1201 

Let /J,„ be a sequence of random probability measures on R converging weakly to a prob- 
ability measure /.i. Let F„ denote the (random) distribution function of for A £ M, 

Fn[\) = ^„(-oo, A]. 

Similarly, let F denote the distribution function of /i. It is not very hard to prove that, for 
all A € E such that /J.{A} = 0, F„(A) converges in distribution to F{\). (We prove this in 
the appendix of the paper, but it is certainly written somewhere in the literature.) 

Let us now apply this general result to our context. The probability measure /i„ is given 
by with c = >/2. It is known to converge to the random measure /Xise- Assume for 
the moment that this measure does not assign a positive weight to any point. Then, with 
the above notation, F„(A) converges in distribution to ^'(A), for all A G R. But, given the 
definition (QJ of 

F„(A) = 1 - Mn(A, oo) = 1 - -^X+{\n^'y^/2) + -^X^{\n^/^ / ^/2), 

n + 1 n + 1 

where the definition of X„ is extended to all reals by X„(x) = if a; e R\Z. By TheoremsEI 
and El the right-hand side converges in distribution to 1 — Y^{\/y/2). Consequently, the 
tail distribution function of the ISE (that is, ^ise(A, oo)) has the same law as F+(A/\/2). 

It remains to prove that /iise does not weight points positively (almost surely). Let A e R. 
Then 

P(AiiSE{A} > 0) = ^ E(m,se{A}) = 0. (50) 
Let e > 0, and let h\,^ be the function plotted on the right-hand side of Figure El Then 

Mise{A} = lim / /iA,e(y)rfA*isE(y) 

= lim / hxe{r{t))dt (by gU) 
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Taking expectations, we obtain 

E(Ai,sE{A}) = IE (^^ tx=r{t)dt^ = ^ P(A = r{t))dt. 

But P(A = r{t)) = for all t e (0, 1) and A, since r{t) has a density with respect to the 
Lebesgue measure for all t. By Ij5()|l . we conclude that /iise does not weight points positively. 

The last statement of Proposition 1201 is then easily proven, using a the symmetry of /Xise 
and the fact that it does not assign a positive probability to any point. 



6. Other tree models and universality 

6.1. Trees with increments 0,±1 

We consider in this section a slight variation on the previous family of trees: the increments 
of the labels along edges may now be 0, ±1. This family of trees has attracted a lot of interest 
in relation to planar maps |6l ISl lTTl [2i| . 

6.1.1. Enumerative results. As above, let Tj = Tj{t) be the generating function of labelled 
trees in which all labels are at most j, counted by their number of edges. Let Sj = Sj{t,u) 
be the generating function of labelled trees, counted by the number of edges (variable t) and 
the number of nodes labelled j (variable u) . Finally, let Rj = Rj (t, u) be the generating 
function of labelled trees, counted by the number of edges and the number of nodes having 
label j at least. As above, it is easy to write an infinite system of equations defining any of 
the famihes Tj , Sj or Rj . The only difference with our first family of trees is that a third case 
arises in the decomposition of trees illustrated by Figure ^ the leftmost child of the root 
may have label j. In particular, the generating function T = T{t) counting plane labelled 
trees now satisfies 

T = 1 + 3tr^ 

while for j > 0, 

= i + i(r,_i + T, + r,+i)T,. (51) 

The equations of Lemmas and El are modified in a similar way. The three infinite systems 
of equations thus obtained can be solved using the same techniques as in Section [2 The 
solutions are expressed in terms of the above series T = T{t) and of the unique formal power 
series Z = Z{t), with constant term 0, satisfying 

(1 + 4Z + Z2)2 ^ ^ 

Observe that T and Z are related by: 

^ _ 1 + 4^ + ^2 



1 + Z + Z^ ' 

We state without proof the counterparts of Propositions [2 El and [3 

Proposition 22 (Trees with small labels j^). Let Tj = Tj(t) he the generating 
function of trees having no label greater than j. Then Tj is algebraic of degree (at most) 2. 
In particular, 

To = 1 - 16t+ IStTo - 27t^To'^. 

Moreover, for all j > —1, 

(l-ZJ+i)(l-ZJ'+4) 



^■(l-ZJ+2)(l_^j+3)^ 



where Z — Z(t) is given by H52|l . 



38 



MIREILLE BOUSQUET-MELOU 



Remarks 

1. As observed in [HI p. 645], there is an "invariant" function attached to equations of the 
form Ij51|l : for j > 0, 

/(T,_i,T,) = /(T,,T,+i) 

where / is now given by 

I{x, y) = xy{l - t{;x + y))-x-y. 

As explained in the remark that follows Propositions 0] and |3 this can be used to derive 
rapidly from lf5T|l the value of Tq. 

2. As was the case for trees with increments ±1, the trees counted by Tq (equivalently, the 
trees having only non-negative labels) are closely related to planar maps. More precisely, 
there is a one-to-one correspondence between non-negative trees having n edges and planar 
maps having n edges (HJ^]. The coefficients of To{t) are also remarkably simple: 

(1 - 12t)3/2 -1 + m ^ 2.3" f2n 
To[t) - = 2^ 



54f2 A^fn + l)(n + 2) 

n>0 



n 



A combinatorial explanation for the algebraicity of Tq is given in Ql] . 

Proposition 23 (The number of nodes labelled j). For any j € Z, the generating 
function Sj = Sj{t,u) that counts labelled trees by the number of edges and the number of 
nodes labelled j is algebraic of degree at most 4 over Q{T,u) (and hence has degree at most 
8 over Q{t, u) ). More precisely, 

,rr c^2 = - ^TiT - 1)(2T + l)So + (T - i fS^^ 
and all the Sj belong to Q{t,u, Sq). Moreover, for all j > 0, 

where Z = Z(t) is given by l|52|l and fi = ^J.{t,u) is the unique formal power series in t 
satisfying 

(1 + z + z^){i + tizfji + ^^z^r 

^ ' (l + Z)2(l-Z)3(l-^2^3) ■ 

The series ^i{t,u) has polynomial coefficients in u, and satisfies 1) = 0. It has degree 4 
over^^{Z,u) and 16 overQ{t,u). 



Proposition 24 (The number of nodes labelled j or more). Let j e Z. The 

function Rj {t, u) = Rj that counts labelled trees by the number of edges and the number of 
nodes labelled j or more is algebraic over Q{t,u), of degree at most 8. It has degree at most 
2 over Q{T,T), where T = T{t) and T = T{tu). More precisely, it belongs to the extension 
of Q(T, T) generated by 



/4(T + T)2 - TT(4 + STT). 
Moreover, for all j > 0, 

{l + iyZ^){l + i^Z3+^) 



R, = T ■ 



^ (l + i/ZJ+i)(l+i/ZJ+2)' 

where Z = Z{t) is given by Ij52|l and v = v(t, u) is a formal power series in t, with polynomial 
coefficients in u, which is algebraic of degree 4 over <Q(u, Z), and of degree 16 over Q{t, u). 
This series satisfies v{t, 1) = 0. The first terms in its expansion are: 

v{t, u) = {u-l)[l + 2,ut+ (15 u + 14 u^) t^ + (104 u + 117 + 83 u^) t^ + 0{t^)^ . 
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6.1.2. Limit laws. We now endow the set of labelled trees having n edges with the uniform 
distribution, and consider the same random variables as for our first family of trees: Af„, 
the largest label, X„(j), the number of nodes having label j, and finally the number 

of nodes having label j at least. 

Again, we can prove that M„n~^/'* converges in law to A^ise/VS, where A^ise is the 
supremum of the support of the ISE, and that for all A S K, the sequence X^{\n^^^)/n 
converges in law to 5ise('\/3A) where gisE is the tail distribution function of the ISE. The 
arguments are the same as for our first class of trees ('Sections K-i.5l and l5.4ll . 

Hence we could just as well have started from the enumerative results of Section Ifi.l.ll 
rather than from those of Section [2l to characterize the laws of iVisE and 5ise(A) (Propo- 
sitions El and EOI) • More remarkably, we have performed on X„(j) an analysis similar to 
that of Sectional and obtained the same local limit law. In other words, for all A > 0, the 
sequence X„([An^/^J)n~'^/^ converges in law to VS/iseIv^A) where /ise is the conjectured 
density of the ISE, given in Conjecture 1151 

In all three cases, the convergence of the moments holds as well. 

6.2. Naturally embedded binary trees 

We study in the section the incomplete binary trees^ carrying their natural labeUing, as 
shown on the right of Figure Such trees are either empty, or have a root, to which a 
left and right subtree (both possibly empty) are attached. A (minor) difference with the 
two previous families of trees is that the main enumeration parameter is now the number of 
nodes rather than the number of edges. 

6.2.1. Enumerative results. Let Tj = Tj{t) be the generating function of (naturally la- 
belled) binary trees in which all labels are at most j, counted by their number of nodes. Let 
5*^ = Sj(t,u) be the generating function of binary trees, counted by the number of nodes 
(variable t) and the number of nodes labelled j (variable u). Finally, let Rj = Rj{t,u) be 
the generating function of binary trees, counted by the number of nodes and the number of 
nodes having label j at least. It is easy to write an infinite system of equations defining any 
of the families Tj, Sj or Rj. The decomposition of trees that was crucial in Section 12 is now 
replaced by the decomposition of Figure The generating function T = T{t) counting 
naturally labelled binary trees satisfies 

r = i + tT^, 

(as it should!) while for j > 0, 

= l + fT,_iT_,-+i. (53) 
Note that the initial condition is now T-i — 1 (accounting for the empty tree). 




Figure 11. The decomposition of naturally labelled binary trees rooted at j. 



'The author has obtained similar, but slightly heavier results for embedded complete binary trees. 
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The equations of Lemmas Eland El respectively become: 
while 

Rj = l + tRj^iRj+i for i > 1, (55) 

and 

R-j {t, u) = Rj+i{tu, 1 /u) for all j £ 7L. (56) 

The three infinite systems of equations thus obtained can be solved using the same tech- 
niques as in Section [2 The solutions are expressed in terms of the above series T = r(t) 
and of the unique formal power series Z = Z{t), with constant term 0, satisfying 

(57) 

Observe that T and Z are related by: 

1 + Z^ 

T = 



1 - Z + 2'2 ■ 

We state without proof the counterparts of Propositions |2l El and Once again, the results 
below are dying for combinatorial explanations! 

Proposition 25. Let Tj = Tj{t) be the generating function of binary trees having no label 
greater than j. Then Tj is algebraic of degree (at most) 2. In particular, 

(1 - 4t)3/2 _ 1 + 8t - 
° " 2t{l+t) ■ 

Moreover, for all j > —1, 

_ (l-Z^+')il-Z^+') 

^ (l_ZJ+4)(l_^J+5)' 

where Z = Z{t) is given by l|57|l . 

It is easy to check that the above series Tj satisfy the equations IjfiHIl and the initial condition 
T_i = 1. The method we used to discover this product form is again borrowed from |^. 

Remark. For this family of trees as well, we have found an "invariant" function attached to 
equations of the form ll53)l : for j > 0, 

/(T,_i,T,-) = /(T„T,+i) 

where 

Tl \ ( ^ , {x^-x-y + y^)t + x + y 

I{x,y)^{x + y)t +^ '—-\ . 

xy xy 

This can be used to derive rapidly from Ij5^i|l the value of Tq . 

Proposition 26 (The number of nodes labelled j). For any j G Z, the generating 
function Sj = Sj (t, u) that counts binary trees by the number of nodes and the number of 
nodes labelled j is algebraic of degree at most 4 over Q(T, u) (and thus has degree at most 8 
over Q{t, u) ). More precisely, 

T\u-i)^ _ (T - i)'^si - 2r5'o(r ^ 1)2(3 - 9r + rr^) + T^jT"^ + r - 1)^ 

u{T-Soy^ ^ (T-l)(5o-l)(r2 + r5o- 5*0)2 

and all the series Sj belong to Q(t,, u, Sq). Moreover, for all j > 0, 

^ {l + fiZ^){l+fiZ^+^) 
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where Z = Z{t) is given by ljR7|l and ji — ^J.{t,u) is the unique formal power series in t 
satisfying 

^ Z{l + ^iZ)^{l + fiZ^)il + fiZ^) 

^ ^(l + Z)2(l + Z + Z2)(l-Z)3(l-^2^5)' 

The series ^J,{t,u) has polynomial coefficients in u, and satisfies ^{t, 1) = 0. It has degree 4 
overQ{Z,u) and 16 over Q(t,,u). 

Comment on the proof. The proof is similar to the proof of Proposition until one 
computes the equation satisfied by Sq. But then, the relation S'o = 1 + tSf does not allow 
us to conclude that 5*1 belongs to Q{t,u,So)- Instead, we compute the algebraic equation 
satisfied by Si. It is found to have degree 4 over Q{u,T). The above relation between Sq 
and Si shows that 5*0 belongs to the extension of Q{t,u) generated by Si. Comparing the 
degrees implies finally that Q{t,u,So) = Q(t,u, Si). Then Ij54|l shows, by induction on j, 
that all the series Sj belong to this field. 

Proposition 27 (The number of nodes labelled j or more). Let j G Z. The generating 
function Rj{t,u) = Rj that counts binary trees by the number of nodes and the number 
of nodes labelled j or more is algebraic over Q{t,u). More precisely, Rq has degree 16 
over Q{t,u) and degree 4 over Q{T,T), with T = T{tu), and all the series Rj belong to 
Q{T,f,Ro) = Q(t,M,i?o). Moreover, for all j > 0, 

_ {l + i.Z^){\ + vZ'^+^) 

where Z = Z{t) is given by Ij52|l and v = v(t, u) is a formal power series in t, with polynomial 
coefficients in u, which is algebraic of degree 8 over Q{u, Z) and 32 over Q{t, u). This series 
satisfies vit, 1) = 0. The first terms in its expansion are: 

iy{t, u) = {u -l)(t+{u + l)t'^ + {2u^ +3u + 3)t^ + 0(f*)) . 

Comment on the proof. The proof is similar to the proof of Proposition [3 until one 
computes the equation satisfied by Rq. One finds that Rq has degree 4 over Q(T, T), and 
degree 16 over Q{t,u). Using Ij56ll . one then derives an equation satisfied by i?i. Strangely 
enough, it turns out that the minimal polynomials of Rq and Ri over Q(T, T) (or over 
Q{t, u)) are the same. The two series are of course different: 

i?o(i,u) = l + tu + u{l + u)t'^ +u{2u'^ + 2u + l)t^ + u{l + u){Au^ +u + 2)t'^ + 0{t'^), 
Riit,u) = l + t + {l + u)t^ + {u^ + 2u + 2)t^ + {l+u){2u^ +u + A)t^ + 0{t^). 

Let P{x) be the minimal polynomial of Rq and Ri over K = Q(T, T). We want to prove 
that Ri belongs to the extension of K generated by Rq. Note that this property does not 
simply follow from the fact that Rq and i?i are conjugate roots of P. For instance, for a 
generic polynomial P of degree 4 over K, with Galois group S4, the four extensions of K 
generated by the roots of P are different. We are going to determine the Galois group of 
our polynomial P, using the general strategy described in (2 p. 141-142]. The resolvent 
cubic of P, which we denote R below, is found to factor into a linear term and a quadratic 
one. Hence the Galois group of R over K has order 2. This impHes that the Galois group 
G of P over K is either the cyclic group of order 4 or the dihedral group of order 8. In the 
former case, the four extensions of K generated by the roots of P coincide (and are equal 
to the splitting field of P) and we are done. In the latter case, there exists a labelling of 
the four roots of P, say Xq, Xi, X2, X^, such that the group G, seen as a subgroup of the 
permutations of {0, 1, 2, 3}, is 

G = {id, (0, 1), (2, 3), (0, 1)(2, 3), (0, 2)(1, 3), (0, 3)(1, 2), (0, 2, 1, 3), (0, 3, 1, 2)}. 

Then 
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• the simple extensions of K generated by the Xi satisfy K(Xo) = K(Xi) and K(X2) = 

K(X3), 

• the root of the resolvent that belongs to K is F = XqXi + 

The root of R that belongs to IK = Q(r, f) is found to be r = t-^/{u{l + u)) + 0{t-'^). 
We already know two roots of P, namely Rq and which are equal to 1 + 0{i). The other 
two roots are respectively of the form Q2 = —t~'^/u + 0{t~^) and Q3 = —t~^/{l + u) + 0{l). 
Prom the value of Y, we conclude that the above properties hold with Xq = Rq, Xi = Ri, 
X2 = Qi and = Q3. In particular, i?o and R\ belong to the same extension of degree 4 

of Q(r, f). 

Then, an induction on j, based on Ij55|l . implies that for all j > 2, the series Rj belongs to 
the extension of IK = Q(T, T) generated by Rq. Since Ri{t,u) = Ro{tu,l/u) and K(i?o) = 
K(i?i), the field K(i?o) is invariant under the transformation A{t,u) 1-^ A{tu,l/u). This 
property, combined with lf5fijl . impHes that for j > 1, the series R-j belongs to K(i?o)- 

6.2.2. Limit laws. We now endow the set of binary trees having n nodes with the uniform 
distribution, and consider the same random variables as for above: M„, the largest label, 
Xn{j), the number of nodes having label j, and X^{j), the number of nodes having label j 
at least. 

Again, we can prove that Mnn"^/"* converges in law to A''ise, where iVisE is the supremum 
of the support of the ISE, and that for all A € M, the sequence X^{Xn^^^)/n converges in 
law to (?ise(A) where gisE is the tail distribution function of the ISE. The arguments are the 
same as for our first class of trees f Sections K-i.51 and I5.4|l . Hence we could just as well have 
started from the enumerative results of Section 16. 2. 11 rather than from those of Section |2l 
to characterize the laws of A^ise and 5ise(A) (Propositions E| and [201) . 

More remarkably, we have performed on Xn{j) an analysis similar to that of Section^ 
and obtained the same local limit law. In other words, for all A > 0, the sequence 
Xn{[\n^^^\)n~^^^ converges in law to /ise(A) where /ise is the conjectured density of the 
ISE, given in Conjecture 1151 

In all three cases, the convergence of the moments holds as well. 
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Appendix: convergence of random distribution functions 

We want to prove the result stated without proof at the beginning of the proof of Section 

Recall that a sequence of real random variables Z„ converges in law to another random 
variable Z if and only if for all a; € R such that F{Z = x) = 0, 

limP(Z„ <x) = F{Z < x). 

n 

This implies the so-called portmanteau inequality: for all a; e K, 

P(Z < x) < liminf P(Z„ < x) < limsupP(Z„ < x) < P{Z < x). (58) 

Let us now use the notation of Section The convergence of /Lt„ to implies that for 
any bounded Lipschitz function / on M [23 p. 71-74]: 

f{x)d^i„{x) / f{x)d^i(x). 
Jr 

Let A e M and let /A,e and gA,e be the functions plotted in Figure El Then 

9x,t{y)dnn{y) < Fn{\) = Mn(-oo, A] < / f x^t{y)d^n{y) ■ 

R ^R 

Hence, for all a; G M, 

fx,,iy)dfin{y) <x] < P(F„(A) <x)<f( f gx,MdMy) < ^ 
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Since /Lt„ converges to fi, and g\^e is a bounded Lipschitz function, 



A similar result holds for the integral involving fx.e- Thus lf58jl implies 

fx.Ay)My) <A< liminf p(^^„(a) < x) < iimsupP(^^„(A) <x)<p[ I gxAy)My) < ^ 



The integral occurring in the rightmost expression of this inequality is bounded from below 
by oo, A — e], while the integral involving /A,e is bounded from above by oo, A + e]. 
Hence 

P(//(-cx), A + e] <x) < liminf P(F„(A) < x) < limsupP(F„(A) < x) < P(/i(-oo, A - e] < x). 

Taking the limit e — > 0+ gives: 

P(Ai(-oo,A] <x) < liminf P(i^„(A) < a;) < limsupP(i^„(A) < a;) < P(/i(-oo,A) < a;). 

K, in addition, the measure fi does not assign a positive probability to A, the rightmost 
expression in the above inequality equals P(^(— oo. A] < a;). The inequality becomes 

P(F(A) <x) < liminf P(F„(A) < x) < limsupP(F„(A) < x) < ¥{F{\) < x), 

so that for all x such that P(F(A) = a;) = 0, 

limP(i^„(A) < a;) = P(i^(A) < x). 

That is, F„(A) converges in law to F{X). 
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